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La presente invention conceme une methode dite de pseudoisation, groupement ou 
« lumping » pour decrire un fluide forme d'un melange de multiples constituants par un 
nombre reduit de constituants (au moins trois), et une methode permettant roperation 
inverse, methode dite de degroupement ou «delumping», c'est-a-dire une methode 
permettant I'obtention de la composition detaillee du fluide a partir de la description par un 
nombre reduit de constituants. 

La presente invention conceme done une methode de modelisation: 

- permettant d'estimer les proprietes ou le comportement des phases liquide et/ou vapeur 
de melanges de multiples composaiits, tels que rencontres in situ dans les reservoirs de 
petrole ou de gaz ou en surface lors de I'exploitation de ces memes gisements, a raide 
d'une representation a nombre reduit de pseudo-composants, 

- et offrant la possibilite de predire en fonction du temps la composition detaillee de 
fluides produits au cours de la production. 

L'ihvention est particulierement utile pour accelerer les calculs lors de la siinulation de 
I'exploitation de gisements souterrains d'hydrocarbures. 

Un tel modele permet a I'ingenieur de gisement de reduire les temps de calcul des 
simulations de comportement de gisement souriiis a une exploitation, tout en conservant 
une bonne modelisation des mecanismes d'echanges entre phases hydrocarbure, et 
d'etablir des profils compositionnels detailles, ceux-ci etant necessaires en particulier au 
dimensionnement previsionnel et a la gestion d'installations de surface, telles que 
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separateurs, unites de traitement, lignes de transport, etc., done utiles en ingenierie de 
surface et de precedes. 

Etat de la technique 

La modelisation des ecoulements dans un reservoir petrolier ou dans un stockage sous- 
terrain est fondee essentiellement sur I'application au reservoir prealablement maille (ou a 
une portion de celui-ci) de la loi de Darcy bien connue decrivant recoulement de fluides 
en milieux poreux, de lois de bilans de matiere dans chaque unite de volume, des relations 
thermodynamiques regissant revolution des proprietes de phase des fluides tels que la 
viscosite, la masse volumique, sur des conditions initiates, sur des conditions aux limites 
de fermeture de la structure, et sur des conditions aux puits. 

Le modele dit « Black Oil », designe ci-apres en abrege par B.O., est Tun des plus utilises 
en simulation petroliere. II permet de decrire un ecoulement tridimensionnel et triphasique 
(eau-huile-gaz) compressible. Les effluents petroliers intervenant dans ce modele sont 
decrits generalement par un constituant eau, et deux constituants pour le fluide de 
gisement, le mot de constituant recouvrant ici la notion de composant (comme H20 pour 
I'eau) et la notion de pseudo-composant (regroupement de composants), Les constituants 
intervenant dans ce modele sont au nombre de trois : un constituant eau (W), un 
constituant hydrocarbure lourd (H), et un constituant hydrocarbure leger (V). Dans un 
modele de type B.O. dit « strict », le constituant (W) n'est present que dans la phase eau, le 
constituant (H) n'est present que dans la phase liquide hydrocarbure (appelee huile ou 
condensat), le constituant (V) se partage entre les phases hydrocarbure liquide et vapeur 
(phase appelee gaz). Un modele B.O. dit « etendu » diffdre d'un modele B.O. « strict » en 
ce que le constituant (H) se partage entre les phases hydrocarbure liquide et vapeur. Mais 
si Fusage des modeles B.O. est applicable a un grand nombre de cas industriels, il n'en 
reste pas moins a deconseiller dans certains cas, en particulier des cas de gisements de gaz 
a condensats soumis a une injection de gaz sec. 

Un autre modele bien connu de simulation dit « compositionnel » est aussi utilise, ou les 
fluides hydrocarbure sont representes par un nombre plus eleve de constituants, au moins 
trois, souvent davantage, le constituant eau n'etant present que dans la phase aqueuse. La 
modelisation des flux de ces fluides mieux detailles se traduit par des temps de calcul 
d'autant plus importants que le nombre de constituants est eleve. 
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Pour que les cailculs de modelisation puissent Stre menes dans un delai raisonnable, les 
fluides en place sont decrits comme constitu^s d'un nombre de composants ou pseudo- 
composants beaucoup plus reduit que le nombre reel de composants. Le passage d'une 
representation detaillee des fluides k une representation a nombre moindre de constituants 
s'effectue par des operations dites de « lumping » (groupement ou regroupement) ou de 
« pseudoi'sation ». Dans ce qui suit, sauf indication contraire, le terme de « pseudoisation » 
sera employe pour toute methode permettant.de reduire le nombre de constituants. 

Differentes m^thodes de pseudoisation ont deja ete proposees pour selectionner, definir les 
pseudo-composants, Tingenieur se trouvant souvent a devoir trouver un compromis entre 
la precision, et le temps (et le cout) des calculs. Pour la simulation de I'exploitation de 
gisements de gaz a condensat soumis a une injection de gaz sec, ce sont des 
representations a quelque 6 ou 8 constituants qui sont en general utilisees, ce qui se traduit 
par des temps de calcul d'autant plus importants qu'il est souvent necessaire de reduire la 
taille des mailles pour limiter I'erreur numerique et par consequent d'augmenter le nombre 
de mailles. Des efforts considerables sont done encore aujourd'hui eonsacres a 
r elaboration de methodes de pseudoisation a I'usage de Tindustrie, qui permettraient de 
reduire encore le nombre de constituants, tout en modelisant le comportement des fluides 
avec precision, et en offrant la ppssibilite de.recouvrer une information compositionnelle 
detaillee. Les operations qui permettent de predire quels auraient ete les resultats. d'une 
simulation de reser\w si Ton avait utilise une modelisation finement detaillee (ou les 
fluides sont representes par un nombre plus eleve de composants) sont connues des gens 
de I'ait sous I'appellation de « delumping », que Ton designera par la suite aussi par 
« eclatement », 

Par le brevet WO 00/37898, on connait une methode de psetidoi'sation, applicable a des 
simulations compositionnelles, basee sur la selection d'un nombre de composants de base 
'dominants' egal au nombre de pseudo-composants que Ton veut obtenir a Tissue de la 
procedure. Dans cette methode, les composants non dominants sont representes dans tous 
les pseudo-composants, et un composant dominant particulier n'est represente que dans un 
seul pseudo-composant. La transformation mathematique a la base de la methode de 
pseudoisation permet, par transformation inverse, de recouvrer Tinformation 
compositionnelle detaillee. La composition d'un pseudo-composant pris en particulier peut 
faire apparaitre des fractions molaires negatives de constituants de base, comme dans 
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Fexemple donne table E7 du brevet en reference. On congoit qu'une telle representation 
perd un certain sens physique des lors que Ton considere isolement un pseudo-composant 
particulier. Cela semble pouvoir poser un probleme de robustesse, lorsque, comme il arrive 
en pratique dans les cas d'injection de gaz, un resultat 'local' de la simulation, par exemple 
5 dans certaines mailles, indique la disparition de Tun ou plusieurs des constituants utilises 
dans la simulation dynamique. Par ailleurs, la mise en oeuvre de I'invention est decrite 
connme necessitant de nombreux calculs iteratifs et un large espace memoire. 

Par des publications anterieures, on connait des methodes de pseudoi'sation, applicables 
egalement a des simulations compositionnelles, oix chaque pseudo-composant est forme en 

10 regroupant plusieurs constituants de base, un constituant de base particulier ne se trouvant 
represente que dans un seul pseudo-composant. Le regroupement peut etre effectue selon 
un choix de criteres fixes a priori, tels que ceux cites dans le brevet mentionne ci-dessus, 
ou par une procedure d'optimisation, par exemple telle que proposee par K. Liu dans le 
papier « Reduce the Number of Components for Compositional Reservoir Simulation », 

15 SPE 66363, presente a la SPE Reservoir Simulation Symposium, Houston, Texas, 11-14 
Fevrier 2001. 

Un papier, bien connu des gens de Tart, ecrit par D.E. Kenyon et G, Alda Behie « Third 
SPE Comparative Solution Project: Gas Cycling of Retrograde Condensate Reservoirs », 
SPE 12278, Journal of Petroleum Technology, Aout 1987, illustre une situation non 

20 exceptionnelle : il arrive de constater, pour un meme cas d'etude, une grande disparite 
dans les resultats de simulations compositionnelles lorsqu'elles sont realisees avec des 
representations de fluides differentes, et qui plus est, avec des logiciels de simulation 
differents. Dans ce papier, le nombre des constituants varie de 5 a 16, et la disparite des 
resultats est en partie imputable aux divers choix de representations compositionnelles des 

25 fluides: la tendance generale observee est que plus le nombre de constituants est faible, 
plus la saturation en liquide hydrocarbure (appele huile ou condensat) dans une maille 
particuliere peut etre sous-estimee, et plus la recuperation d'huile en surface est alors 
surestimee. 

Pour recouvrer I'infoiination compositionnelle detaillee au cours d'une simulation 
30 compositionnelle de reservoir, on peut utiliser des methodes de delumping telles que celle 
exposee dans le brevet WO 99/42937 et aussi dans le papier de C. Leibovici et J. Barker 
« A Method for Delumping the Results of a Compositional Reservoir Simulation » SPE 
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"^ 49068, presente a la SPE Annual Technical Conference and Exhibition New-Orleans, 27- 

, 30 Septembre 1998. La methode permet de prevoir revolution de la composition detaillee 
au cours du temps a partir des calculs, en particulier de coefficients d'equilibre, effectues 
dans une simulation de type compositionnel de fluides decrits par un certain nombre reduit 
5 de pseudo-composants, le nombre de composants etant au moins de trois. 

Une representation de type Black-Oil peut etre consideree comme issue d'une operation de 
pseudoi'sation particuliere foumissant deux pseudo-composants. La composition detaillee 
de chacun des pseudo-composants, qui n'est pas utile a la construction de la 
representation, n*est pas connue a priori, sans que ceci soit un obstacle redhibitoire au 
10 recouvrement d'une information compositionnelle detaillee, par une operation de 
delumping. Ainsi, par le brevet FR 00/09008, on connaat. une methode qui permet de 
prevoir revolution de la composition detaillee au cours du temps a partir des calculs 
effectues dans une simulation dynamique de type Black-Oil. 

Le principe de Tetape d'eclatement (ou delumping) des brevets WO 99/42937^ et FR 
15 00/09008 est de determiner, - a partir des calculs effectues pendant la simulation, avec la 
representation thermodynamique regroupee (compositionnelle ou BO), dans chaque maille 
et a chaque pas de temps -, le coefficient ADo et les n coefficients ADp (soit -n+1 
coefficients, n etant le nombre de parametres de Tequation d'6tat) d'une equation generale 
connue, publiee anterieurement dans un papier de C.F. Leibovici, E.H. Stenby, K. 
. 20 Knudsen, « A Consistent Procedure for Pseudo-Component Delumping », Fluid Phase 
Equilibria, 1996, 117, 225-232, : 

(1) Ln(Ki) = ADo + X ADp Opi 

ou Ki est la constante d'equilibre du constituant i et les Ilpi sont des parametres fixes de 
caracterisation du constituant i dans Tequation d'etat pour une representation 
25 thermodynamique donnee. 

Une fois le coefficient ADo et les n coefficients ADp determines, on^les utilise pour calculer 
les constantes d'equilibre des constituants de la representation thermodynamique detaillee 
(Nrb composants) en appHquant T equation (1) aux Nrb composants avec leurs propres 
parametres fixes de caracterisation dans la representation thermodynamique detaillee. 
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L'un des interets de cette methode est que, dans I'etape d'eclatement (« delumping »), il 
n'est pas necessaire de resoudre les Nrb equations d'^quilibre associees a I'equation d'etat 
(equations qui traduisent I'egalite des fugacites de chaque constituant dans chaque phase) 
aux differents pas de temps de la simulation dynamique des ecoulements, ce qiri constitue 
une economie de temps de calcul. 

Le papier de W.H. Goldthorpe « Simulation of Gas Injection Processes in Gas-Condensate 
Reservoirs Using a Binary Pseudo-Component Representation », SPE 19470, present^ a la 
SPE Asia-Pacific Conference, Sydney, Australie, 13-15 Septembre 1989, illustre les 
resultats de simulation que Ton peut obtenir avec une modelisation Black-Oil avancee 
dans un cas d' exploitation d'un reservoir de gaz k condensat par un precede d'injection de 
gaz. Le cas pris comme exemple est issu de la publication ci-dessus mentionn6e de D.E. 
Kenyon et G. Alda Behie. Au vu de la disparity des resultats dans cette publication, les 
resultats de W.H. Goldthorpe avec une representation Black-Oil, jugee par comparaison 
avec les resultats d'une simulation effectu^e avec une representation detaillee, apparaissent 
bien plus coh6rents, mais on peut constater dans la figure 8 du papier que la solution n'est 
pas satisfaisante dans la phase de re vaporisation car la saturation en huile dans une maille 
particuliere (la meme que dans le papier de reference) est tres differente de celle de la 
prediction compositionnelle detaillee, et semble par ailleurs etre tronqude de valeurs 
negatives pendant huit annees d' exploitation simulee. 

L'etat de la technique antdrieure est done que I'on ne dispose pas de methode de 
pseudoifsation a la fois simple et robuste: 

- permettant de reduire k trois le nombre de pseudo-composants dans des simulations 
compositionnelles, pour beneficier de temps de calcul notablement reduits, notamment 
pour les simulations de cas d'injection de gaz avec effets de re vaporisation, 
difficilement traitables en Black-Oil, 

- garantissant un sens physique aux resultats des simulations compositionnelles et, par 
voie de consequence, aux resultats d' operations associees de delumping. 

La methode selon l'invention 

La methode de pseudoisation selon l'invention permet d'estimer les proprietes ou le ■ 
comportement des phases hydrocarbure liquide et/ou vapeur, a partir de donnees relatives 
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^ un ensemble de reference constitue de melanges hydrocarbones dans une serie d'etats 
thermodynamiques decoulant de conditions d* exploitation rencontrees, ou prevues comme 
telles, pour des gisements souterrains d'hydrocarbures. Elle comporte les etapes suivantes : 

on regroupe chacun des dits melanges hydrocarbones en au moins trois constituants 
(V, I, H), aucun de ces constituants ne correspondant a une selection particuli^re de 
composants ou pseudo-compos ants de base qui serviraient k une description 
compositionnelle detaillee des fluides, en considerant que les phases gaz issues de la 
separation en conditions de surface de chacun des melanges hydrocarbones sent des 
melanges dont le troisieme constituant (H) est exclu, et que les phases huile issues de 
la separation en conditions de surface de chacun des melanges hydrocarbones sont des 
melanges dont le premier constituant (V) est exclu ; 

- on determine par bilan de matiere les compositions des produits de separation 
comprenant, pour les produits gazeux, au moins le premier et le deuxieme constituant , 
(V, I) en proportions variables, et pour les produits Hquides, au moins le deuxieme et le 
troisieme constituant (I, H) en proportions variables ; et 

on determine la composition a au moins trois constituants de chaque melange 
hydrocarbone de T ensemble de reference par combinaison des produits de sa 
separation au prorata des quantites de chaque produit de separation. ^' 

Suivant un mode de mise en oeuvre, on regroupe chacun des melanges hydrocarbones en 
seulement trois constituants (V, I, H), les phases gaz issues de la dite separation sont des 
melanges en proportions variables du premier constituant (V) et du deuxieme constituant 
(I), on regroupe les phases huile issues de la dite separation sont . des melanges en 
proportions variables du deuxieme constituant (I) et du troisieme constituant (H), et on 
determine la cornposition a trois constituants. 

Les conditions de surface sont celles rencontrees ou prevues lors de 1' exploitation du 
gisement mais elles sont redefinissables selon le contexte de Tetude menee, et peuvent 
done etre differentes des conditions de surface rencontrees ou prevues lors de 
Texploitation du gisement. j . 
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Suivant un detail de mdse en oeuvre, le bilan de matiere est un bilan massique et on attribue 
a chacun des trois constituants (V, I, H) une masse molaire apres une analyse quantitative 
des masses mol aires des produits de separation de Tensemble de reference. 

Suivant une variante de mise en ceuvre, on definit les donnees necessaires aux calculs 
d'equilibre et a la modelisation des proprietes de phases dans la representation groupee, en 
utilisant les compositions des phases dans la representation groupee et des donnees 
connues, ou estimees, a priori portant au moins sur la densite et la viscosite de phases huile 
et gaz en equilibre appartenant a I'ensemble de reference. 

Lorsque les dites donnees incluent des donnees compositionnelles detaillees des phases 
prealablement representees par une equation d'etat « detaillee », on definit les parametres 
d'une premiere (le cas echeant unique) equation d'etat de la representation groupee, utile a 
la modelisation des proprietes de phases, en utilisant ces donnees compositionnelles. 

Suivant une autre variante preferee de mise en oeuvre, on ajuste les parametres d'une 
seconde equation d'etat de la representation groupee, utile aux calculs d'equihbre, pour 
reproduire les coefficients d'equilibre de la representation groupee. 

Pour r^aliser cet ajustement, on utilise par exemple les parametres par constituant de la 
representation groupee dans une equation d'etat utile aux calculs des proprietes des phases. 

Suivant une autre variante de mise en oeuvre utile en particulier k des calculs ulterieurs 
d'eclatement ou « delumping », on determine les coefficients d'equilibre des fluides dans 
une representation compositionnelle detaillee, a partir de variables et/ou parametres 
intervenant dans le calcul des proprietes de phases, des lors que les parametres utiles aux 
calculs des proprietes de phases dans la representation groupee ont etc estimes de fagon a 
reproduire les parametres des phases dans I'equation d'etat de la description 
compositionnelle detaillee. 

La methode peut egalement comporter des etapes d'eclatement ou delumping pour predire 
en fonction du temps et dans au moins une zone thermodynamique, une composition 
detaillee d'un fluide contenu dans un gisement d'hydrocarbures ou produit par au moins 
un puits, ces etapes etant par exemple les suivantes : 
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on represente le reservoir sous la forme d'un reseau de mailles (m) dont chacune 
constitue un volume elementaire rempli de fluides sous forme de une ou plusieurs 
phases, avec au moins une phase non aqueuse ; 

on d^finit, pour chaque zone bu domaine thermodynamique, les fluides par une 
representation detaillee de base, de fagon a determiner la quantity de chaque 
constituant de base (i) dans chaque phase hydrocarbure dans chaque maille (m) k 
. r instant defini comme initial pour le calcul de delumping ; 

par zone thermodynamique pour laquelle on adopte une representation groupee des 
fluides, on determine une equation d'etat construite, prealablement a la simulation 
dynamique de reser\^oir avec la representation groupee, pour reproduire les parametres 
de phases, dans I'equation d'etat de la representation detaillee, des fluides 
hydrocarbones au cours de chemins thermodynamiques tenus pour representatifs de 
ceux que vont suivre les fluides hydrocarbones pendant la simulation maillee ; 

on realise a un pas de temps t une simulation compositionnelle avec un nombre Hmite 
de constituants ou les proprietes des phases sont calculees par une equationnd'etat, 
ladite simulation permettant de calculer au moins dans chaque maille (m) et a des pas 
de temps consecutifs, une pression pour une phase hydrocarbure, la temperature 
lorsque celle-ci varie, les debits des phases entre mailles, entre mailles perforees et 
puits, et les valeurs de parametres et/ou proprietes de phase intervenarit dans 
. r expression formelle des coefficients d'equilibre de la representation detaillee, et on 
memorise ces di verses quantit^s ; 

on estime au pas de temps suivant (t+1) la fraction molaire de chaque constituant i 
dans la composition detaillee globale du fluide hydrocarbone dans la maille (m) par 
bilan matiere sur la maille (m) ; 

on determine en utilisant les quantites memorisees, au meme pas de temps (t+l) et 
dans chaque maille (m), les coefficients d'equilibre de chaque constituant (i) dans la 
representation detaillee ; 

on determine, au meme pas de temps (t+l), la fraction vaponsee dans chaque maille 
(m) ; et 
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- on estime la composition detaillee de chaque phase hydrocarbonee, au meme pas de 
temps (t-f 1) et dans chaque maille (m) en particulier sous la forme des nombres de 
moles de chaque constituant i dans la representation detaillee dans chaque phase au 
meme pas de temps (t+1). 

En d'autres termes, on peut aussi defmir Fessentiel de la methode de pseudoi'sation comme 
suit : 

- Ayant a disposition des donnees, au moins de densites et viscositds, concemant des 
phases huile et gaz en equilibre «in situ » pour un certain nombre d'etats 
thermodynamiques 'e' dans des conditions de pression et temperature et T^ 

- connaissant directement ou par Tintermediaire de calculs preliminaires, les masses 
molaires des phases *huile' en provenant de la separation de Fhuile in situ, de la 
separation du gaz in situ, notees respectivement (MMoo) et (MMgo), les masses 
molaires des phases gaz en surface provenant de la separation de I'huile in situ, de la 
separation du gaz in situ, notees respectivement (MMog) et (MMgg), et les fractions 
molaires vapeur (ou fraction molaire de phase gaz). Boo pour la separation de I'huile in 

e 

Situ, Bgo pour la separation du gaz in situ, 

- on determine les compositions des phases « en surface », et par voie de consequence 
les coefficients d'equilibre, dans la representation groupee utilisant les constituants (V, 
I, H) avec Thypothese que les phases gaz en surface peuvent contenir les constituants 
(V) et (I), a Fexclusion du constituant (H), et que les phases huile en surface peuvent 
contenir les constituants (I) et (H), a Fexclusion du constituant (V). 

- on determine les compositions des phases « in situ » dans la representation groupee 
utilisant les constituants (V, I, H) par recombinaison des compositions des phases en 
surface en utilisant les fractions molaires vapeur 0o6, 9go, ces determinations 
foumissant les coefficients d'equilibre des constituants (V, I, H) dans les conditions de 
reservoir « in-situ », 

- on utilise les donnees connues a priori et les compositions des phases dans la 
representation groupee pour definir les entrees necessaires a la modelisation des 
proprietes de phases et aux calculs d^equilibre dans la representation gi*oupee et, plus 
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particulierement lorsque les donnees a priori incluent des donnees compositionnelles 
detaillees des phases prealablement representees par une equation d'etat, on utilise ces 
informations additionnelles pour definir les parametres d'une equation d'etat de la 
representation groupee utile en premier lieu a la modelisation des proprietes de phases, 
5 le cas echeant pour definir les parametres d'une seconde equation d'etat de la 

representation groupee utile aux calculs d'equilibre. 

Les constituants V, I, H s'apparentent respectivement a un constituant volatile, un 
constituant intermediaire et un constituant lourd, sans qu'aucun d'entre eux ne corresponde 
a une selection particuliere de composants ou pseudorcomposants de base qui serviraient a 
10 une description compositionnelle detaillee des fluides. 

L'idee essentielle mise en oeuvre pour le degroupage ou delumping des simulations 
compositionnelles est de calculer les coefficients d'equilibre des fluides dans la 
representation compositionnelle detaillee a partir de variables et/ou parametres intervenant- 
dans le calcul des proprietes de phases des lors que les parametres utiles aux calculs des 
15 proprietes de phases dans la representation groupee ont 6t6 estimes de fa§on a reproduire 
les parametres des phases dans I'equation d'etat de la description compositionnelle 
detaill6e 

Presentation des figures ^ ^ 

D'autres caracteristiques et avantages de la methode selon I'invention, apparaltront a la 
20 lecture de la description ci-apres d'un exemple non limitatif de realisation, en se referant 
aux dessins annexes presentes ou : 

la figure 1 schematise la separation distincte d*une phase huile et d'une phase gaz en 
equihbre dans le reservoir ; 

- les figures 2 a 5 comparent les resultats d'une simulation compositionnelle avec une 
25 representation a trois pseudo-composants (resultats d^signes par 'temaire'), obtenue 

par application de la methode selon I'invention, a ceux d'une simulation 
compositionnelle detaillee a seize constituants de base (resultats desisnes par 
'reference ; 

- les figures 6 et 7 montrent, pour les representations detaillee et temaire, respectivement 
30. les variations du terme d'attraction et du co- volume (exprimes sous forme 
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adimensionnelle) des phases huile et gaz en fonction de la pression au cours de la 
simulation thermodynamique d'une operation de depressurisation a volunae constant a 
la temperature de reservoir ; 

- les figures 8.1 a 8.16 d'une part, et les figures 9.1 a 9.16 d' autre part, montrent, par 
5 constituant, respectivement la composition detaillee du gaz et de I'huile, au cours de la 

meme operation de depressurisation, telles qu'obtenues a partir de la representation de 
reference, et a partir de I'eclatement de la representation a trois pseudo-composants ; 

la figure 10 recapitule les resultats illustres aux figures 8.1 a 8.16 et 9.1 a 9.16, en 
montrant I'erreur absolue maximale obtenue sur la teneur en chaque constituant durant 
10 la depressurisation de 237.4 a 69 bars en fonction de la teneur du constituant dans la 

phase a la pression de rosee du fluide initial ; 

les figures 11.1 a 11.16 d'une part, et les figures 12.1 a 12.16 d'autre part, montrent, 
par constituant, respectivement la composition detaillee du gaz et de Thuile, au cours 
de la simulation thermodynamique d'une operation de vaporisation d'un condensat par 
15 injection d'un gaz sec a la pression constante de 169 bars, telles qu'obtenues a partir de. 

la representation de reference, et a partir de I'eclatement de la representation a trois 
pseudo-composants ; et 

la figures 13 recapitule les resultats illustres aux figures 11.1 a 11.16 et 12.1 a 12.16, 
en montrant I'erreur absolue maximale obtenue sur la teneur en chaque constituant 
20 pour I'operation de vaporisation en fonction de la teneur du constituant dans la phase 

au demarrage de I'operation de vaporisation. 

Description detaillee 

Preambule 

La modelisation du comportement de melanges fluides necessite un mode de definition des 
25 constituants utilises pour decrire la composition des fluides et comporte : 

des calculs d'equilibre thermodynamique, pour savoir si, dans des conditions donnees 
de pression et de tempei*ature, I'etat d'equilibre d'un melange, caracterisd par une 
composition globale, est monophasique, ou si le melange se separe en differentes 
phases, le calcul d'equilibre permettant alors de connaitre les proportions de chaque 
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phase, et la composition de chaque phase. Un calcul d'^quilibre sur un melange, dont 
la composition globale est decrite par les fractions molaires z\ de chaque constituant 
(i), conduit a la determination de la fraction vapeur 9 et a la determination de la 
composition de chaque phase. La fraction vapeur 9 est determinee en general en 
5 resolvant Tequation de Rachford-Rice connue des specialistes qui est appliquee par 

exemple dans le document suivant : 

- Rachford, H.HJr et Rice, J.D.; « Procedure for Use of Electronic Digital Computers in 
Calculating Hash Vaporization Hydrocarbon Equilibrium", J. Pet. Technol., 1952, 14, 
19, 

les compositions des phases 'huile' et *gaz*, decrites respectivement par les fractions 
molaires Xj et yi, etant determinees par: 

5i Kj Zj 

Kp) ^i-i+(Ki~l)9 y*"l+(Ki-l)9 

Une phase hydrocarbure est dite sous-saturee a une pression et une temperature T donnee 
15 lorsqu'elle ne se trouve pas en .6quilibre avec une autre phase hydrocarbure. Selon une 
convention souvent admise par les gens de I'art, une phase hydrocarbure sous-saturee est 
dite/huile' (ou 'condensat') lorsque sa temperature critique est inferieure a la temperature 
T, et est dite 'gaz' lorsque sa temperature critique est superieure a la temperature T, 
Lorsque sa temperature critique est egale a la temperature T, la phase pent etre qualifiee, 
20 en principe indifferemment, de 'huile' ou 'gaz', et est done consideree comrne rattachee ^ 
Tune ou Tautre des qualites *huile' ou 'gaz'. 

des calculs de proprietes des phases telles que masse volumique, viscosite, 

le cas echeant, des calculs, dits de delumping, permettant la restitution d'une 
information compositionnelle detaillee lorsque les fluides ont ete modelises de fagon 
25 non detaillee. ' 

Les differentes etapes de la methode de modelisation proposee sont presentees ci-apres, la 
presentation etant structuree en cinq sous-chapitres: 
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mode de definition des constituants, 
proprietes des phases, 
equilibres thermodynamiques, 
delumping, 

remarques concemant la g6n6ralisation de rapproche. 



Mode de definition des constituants 

a) On considere un ensemble «E», dit de reference, de fluides non aqueux, dits 
melanges hydrocarbure mSme s'ils peuvent contenir certains composants autres que des 
hydrocarbures tel que I'azote, le dioxyde de carbone, I'anhydride sulfureux, susceptibles 
de correspondre : 

- aux phases fluides hydrocarbure sous-saturees trouvees, supposees ou predites dans 
une ou plusieurs zones thermodynamiques du reservoir consider6 ; 

- aux phases fluides hydrocarbure en equilibre trouvees, supposees ou predites dans une 
ou plusieurs zones thermodynamiques du reservoir considere ; 

- aux phases fluides hydrocarbure correspondant aux fluides hydrocarbure injectes (ou 
supposes I'etre) dans le reservoir. Dans les cas pratiques les plus frequents, les fluides 
d'injection sont des gaz contenant principalement des hydrocarbures 16gers et 
intermddiaires. 

b) Lorsque Ton dispose des donnees n^cessaires, on definit une representation 
detaillee a Nh, composants et/ou pseudo-composants de base pour caracteriser les 
melanges fluides. Cette situation dtant tnes courante, on place la suite de I'expose dans une 
telle situation car cela permet d'introduire progressivement les conventions d'usage utiles 
a la description de la methode. 

Le nombre N^b etant typiquement sup6rieur k 10, chaque constituant est decrit par un 
certain nombre de parametres, dont sa masse molaire MMi, I'ihdice i etant le num^ro 
attribue au constituant ( i=l,..., Hb). Les fractions molaires Xi et y, decrivent 
respectivement la composition des phases hydrocarbure liquide ('huile') et vapeur ('gaz') 
dans la representation de base. 
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Un element 'e' de 1' ensemble « E » est alors decrit par : une pression P , une temperature 
T^ et: 

(2 X Nrb) fractions molaires (xi xna) (yi yNrb) lorsque deux phases hydrocarbure 
en equilibre coexistent, 

e e e e 

5 - ou (Nrb) fractions molaires (xi XNrb) ou (yj ywrb) lorsque une seule phase 
hydrocarbure (huile ou gaz) sous-saturee est pr6sente. 

Dans ce qui suit, on ecrit les equations pour un element 'e' oil deux phases liquide et 
vapeur coexistent. Pour un element *e', ou une seule phase sous-saturee est presente, il 
suffit de selectionner le jeu d'equations relatives a la seule phase representee. 

10 c) On considere que Ton dispose, pour chaque phase de Fensemble «E» de 
reference, de donnees de fluides susceptibles de correspondre aux corripositions des phases 
fluides hydrocarbure trouvees, supposees ou predites correspondant aux melanges fluides 
produits par une chaine de separation, la chaine de separation pouvant etre reduite a une 
seule detente ou comporter un etage ou plusieurs etages intermediaires de separation. Une 

15 chaine de separation est ici, comme il est d'usage dans ringenierie de gisement, 
consideree comme definie par les valeurs de pressions et temperatures des etages de 
separation. H est pratique de considerer la chaine de separation existante. ou prevue par 
laquelle le melange fluide est ou est suppose etre produit. On peut etre amene, dans des cas 
pratiques, a considerer differentes chames de separation, par exemple une chaine de 

20 separation differente par zone thermodynamique, ou un changement de chaine de 
separation dans une zone thermodynamique donnee. Dans ce qui suit, on simplifie d'abord 
la presentation en ne consid6rant qu'une chaine de separation identique pour tous les 
elements de Tensemble « E », puis on considere le cas de plusieurs chaines de separation. 

6 C 

Les produits de separation d'une phase huile de composition (xi x^rb) sont en general une 
25 phase liquide (phase huile) et une phase vapeur (phase gaz), dont les masses molaires sont 
supposees connues et designees respectivement par MMqo et MMqg, dans les proportions 
molaires supposees connues de 9o de phase vapeur et (1 - 8o) def phase liquide. Go peut 
prendre toutes les valeurs reelles entre 0.0 et 1,0, homes incluses. L'indice o (en caractere 
gras) rappellequ'il s'agit des produits de sepai*ation de la phase huile. Les indices o et G 
30 (en caractere non gras) donnent la qualite huile ou gaz du produit de separation. 
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Les produits de separation d'une phase gaz de composition (yi yNrb ) sont en general une 
phase liquide (phase huile) et une phase vapeur (phase gaz), dont les masses molaires sont 

e e 

supposees connues et designees respectivement jpar MMoc MMgg. dans les proportions 

molaires supposees connues de 9g de phase vapeur et (1 - 9g) de phase liquide. Gg pent 
prendre toutes les valeurs reelles entre 0.0 et 1.0, homes incluses. L'indice G (en caractere 
gras) rappelle qu'il s'agit des produits de separation de la phase gaz. 

La Fig. 1 schematise la separation distincte d'une phase huile et d'une phase gaz en 
equilibre dans le reservoir dans la situation generale ou la separation (simple ou complexe) 
d'un melange initialement monophasique, liquide (huile) ou gaz, produit deux phases, 
Tune gazeuse, 1' autre liquide. 

Les relations suivantes traduisent la conservation de la masse lors des operations de 
separation respectivement de I'huile et/ou du gaz de I'element 'e': 

(4) r MmS = (1 - eS) MMoS + Go MMoS 

MMg = (1 - Gg) MMgo + Gg MMgg 

e e 

ou MMo et MMg sont les masses molaires des phases huile et gaz d*un element 'e'. 

d) La methode de pseudoisation proposee consiste a decrire chaque phase 
hydrocarbonee de I' ensemble « E » en au moins trois constituants (V, I, H). 

On definit les pseudo-composants en considerant une hypothese designee ci-apres par 
« M », selon laquelle les phases gaz issues de la separation en conditions de surface de 
chacun des melanges hydrocarbones sont des melanges dont le troisieme constituant (H) 
est exclu, et que les phases huile issues de la separation en conditions de surface de chacun 
des melanges hydrocarbones sont des melanges dont le premier constituant (V) est exclu. 

Les constituants (V), (I), (H) s'apparentent respectivement a un constituant volatile, 
intermediaire et lourd, sans qu'aucun d'entre eux ne corresponde a une selection 
patticuliere de composants ou pseudo-composants de base qui serviraient a une description 
compositionnelle detaillee des fluides. La modelisation associee est designee ci-apres en 
abrege par BRO, pour 'BROWN OIL', pour evocation de la couleur de I'huile qui passe 
du noir au marron avec un contenu croissant en constituants intermediaires. 
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En un premier temps, pour faciliter Texpose, on considere selon une interpretation 
restrictive de Thypothese maitresse «M» que les melanges produits par la chaine de 
separation sont: 

un melange en proportions variables des pseudo constituants volatile (V) et 
5 intermediaire (I) pour un produit gazeux (la masse molaire du constituaiit V etant 

inferieure a celle du constituant I) ; . 

et un melange en proportions variables des pseudo constituants intermediaire (I) et 
lourd (H) pour un produit liquide Qa niasse molaire du constituant I etant inferieure a 
celle du constituant H). 

10 La restriction sur Thypothese « M » sera levee dans le paragraphe tt). 

Les relations suivantes s'appliquent done: 

(5) ( MMoo = ^oi MMi + xoH MMh ; xoi + xqh = 1 

MMoG = yov MMv + yoi MMi ; yov + yoi = 1 

MMgo = xgi MMi + xgh MMh ; xgi + xgh = 1 

MMgg = ycv MMv + ySi MMi ; ySv -I- yci = 1 
Pour conserver un sens physique a la representation temaire, il est necessaire que: 
les masses molaires MMv, MMi, MMh soient definies positives ; 

eeeeeeee. 

les fractions molaires xqi , xqh , yov , yoi , xgi , xgh > ycv , yci soient comprises entre 
15 0.0 et LO, bomes incluses. 

II en decoule que le choix des masses molaires MMv, MMi, MMh doit respecter les 
inegalites suivantes: 



MMv ^ Min {MMgg, MMog} 



l(eeE) 



Max {MMgg, MMog} I ^ MM, < Min {MMgo,MMoo}| 

|(eeE) |(eeE) 



e _ _ ^ e 



20 . ; MMh ^ Max {MMgo, MMoo) 



(eeE) 
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Dans les equations precedentes, les accolades designent un ensemble de valeurs sur 
lesquelles on recherche une valeur minimum (Min) ou maximum (Max). 

Les inegalites etant non strictes, I'absence de Tun des constituants (V,I) dans un produit de 
separation gazeux donne est possible mais ne constitue pas une regie generale et, de meme, 
r absence de Tun des constituants (I,H) dans un produit de separation liquide donne est 
possible mais ne constitue pas une regie generale. 

e) Les valeurs des masses molaires des constituants (V), (I), (H) doivent etre choisies 
en respectant les conditions ci-dessus et etre toutes differentes les unes des autres. Dans 
Texemple foumi a titre d' illustration, on a pose, par analogie avec une demarche Black- 
Oil, et sans que ce soit un exemple limitatif de realisation: 



MMv = Min {MMgg, MMqg} 



(eeE) 



MMh = Max {MMgo. MMqo} | 



(eeE) 



L'exemple d' illustration, largement utilise dans la suite de la description sous la 
designation de cas *SPE3% correspond au premier cas test de la publication deja citee de 
D.E. Kenyon et G. Alda Behie d'exploitation d'un gisement de gaz a condensat par 
15 cyclage de gaz. 

f) Une fois les masses molaires des constituants (V), (I) et (H) definies, la 
composition des fluides issus de la chaine de separation consideree en c) se deduisent, par 
les hypotheses du paragraphe d) et les equations (5), par: 

(6) 



Xov 


= 0 




xoi = 


MMh- 
MMh 


MMoo 
-MMi 


e 

; xoH 


MMoo - 
~ MMh- 


-MMi 
MMi 


yov 


MMi - 
~ MMi 


-MMoG 
-MMv 


; yoi = 


MMoG 
MMi- 


-MMv 
■MMv 


yoH 


= 0 




Xgv 


= 0 






MMh- 


MMgo 


e 


MMgo - 


-MMi 








xgi = 


MMh- 


-MMi 




~ MMh- 


MMi 


yov 


MMi- 
~ MMi 


-MMgg 
-MMv 


; yGi = 


MMgg 
MMi- 


-MMv 
MMv 


■> 


= 0 
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g) Pour obtenir la composition des phases de 1' element 'e' dans la representation a 
trois constituants, on utilise les relations suivantes (dans lesquelles une notation en 
caracteres gras est adoptee pour designer les fractions molaires des fluides in situ) qui 
ti-aduisent la conservation molaire lors des operations de separation respectivement de 
I'huile et/ou du gaz de I'element 'e' : 



(7) 



Xk = (1 - 9o) XOK + ©O yOK 



K = V,I,H 



Yk = (1 - Qg) xgk + 9g Ygk 

En introduisant les relations (6) dans le jeu d'equations (7), on obtient dans le reservoir 
pour une phase huile de Telement *e' la composition suivante: 



(8) 



/■ e 
e MMi - MMqg 

^^"^^ MMi-MMv 

c_ e MMh -MMqo e MMqg MMy 

x,-(l~eo) mMh-MMi "^^^ MMi-MMv 

e e MMqo-MMi 

I XH = (l-eo) MMa-MMi 



et pour une phase gaz: 
(9) ^ 



e _ e MMi -MMgG 

®^ MMi - MMv 

e_ e MMh "MMgo c MMqg - MMy 

yi = (l-0G) ]viMH-MMi "^^^ MMi-MMv 



e e MMgO-MMi 

yH=-(l -eG) MMh -MMi 



h) A Tissue des etapes a) a ,g), on dispose done, pour une utilisation ulterieure, des 
donnees minimales pour definir les constituants (V), (I), (H) a savoir leurs masses 
molaires. / 
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On note que, pour I'ensemble des elements 'e' correspondant a des etats d'equilibre entre 
phases hydrocarbure, on dispose, d'apres les jeux d' equations (8) et (9), des coefficients 
d'equilibre de chacun des constituants (V), (I), (H): 



(10) 



f 



rr^ Iv 9g MMi - MMnn 

IVv = e = ^ e ~r 

xv Bo MMi - MMoG 

_ e MMh-MMoo e MMgg-MMv 
j^e_yi_ ^ "^-^ MMh-MMi '^^^ MMi-MMy 

Xl _ e MMh-MMoA e MMqg-MMv 

MMh - MM, MMi - MMy 



e ^ 1 - Qg MMgo - MMi 

IVh = e = ^ ^ e e 

Xh 1 - Go MMoo - MMi 



On note 6galement que le jeu d'equations (6) definit les coefficients d'equilibre pour les 
conditions de sortie, en pression et temperature, de la chaine de separation consid6ree dans 
le paragraphe c): 



(11) 



r e 

JS^V — e = 

xov 



. _ _ ^^OG ~ ^MMy MMh - MMi 
xoi ~ MMh - MMoo - 



^2:[Sv^^ I^^^ yoi MMgg - MMy MMh - MMi 
~ xgi ~ MMh - MMgo ^"^i ~ MMy 



Key = e 
xgv 



; KoH = = 0 

XOH 



> tMJH — e = U 

xgh 



i) Le jeu d'equations (7), ecrit pour respectivement les constituants (V) et (H), donne 
les relations suivantes: 

xy = Go yov ; xh = (1 - Go) x^h 

yv = Gg ycv ; yn = (1 - Gg) xgh 
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doht il decoule: 



(13) 



e 

e Xv 
yOV = -.e 



e 

e Xh 
XOH = 



l-Qo 



e Yh 
Xgh =■ ^ „ e 



1 -0G 



puis: 



(14) 



e 

e ^ Xv 



Xoi = 1 - 



l-Go 



e , Yh 



1-4 

^ 



1 - 



e 

Xh 



Kgi = ' 



1 _^ 



1 - 



c 



1-e, 



G 



Le jeu d*equations (14) montre que, des lors que les masses molaires des constituahts (V), 
(I), (H) sont fix^es, et que les compositions des phases d'un element 'e' dans la 
representation temaire ont ete etablies, les coefficients d'equilibre du constituant (I) 
peuvent etre obtenus a paitir de la seule connaissance des fractions molaires des 
constituants (V) et (H) dans les phases huile et gaz de Teleraent 'e' et des fractions vapeur 

e e 

6o et 9g en sortie de chaine de separation pour chacune des phases, Une premiere 
impilication est que, en cas de chaines de separation multiples, il est possible de generer les 
donnees de coefficients d'equilibre du constituant (I), relative a une autre chaine de 
separation (ci-apres designee par CS_B) que celle utilisee pour obtenir les compositions 
huile et gaz de Telement 'e* (ci-apres designee par CS_A). Ceci est schematise dans le 
tableau suivant : 

f 
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Melange 
detaille 


Operation 


Donnees utilisees 


Melange 


Coefficients 

rl'pmnlihrf> R 

VJl wkJ tiilJlUI C V»-0 xJ 




CS_A 

> 


{MMoo, MMoG, 9o}cs-A 


{Xv, Xi, Xh}cs_a 


e 

- Xv CS_A 
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KqI = e 

yH CS_A 


CS_B 

> 


{6g}csjb 
> 


> 


1 - 9g cs_b 



Les conditions necessaires a une telle procedure sont le respect des inegalites suivantes: 
yoi cs_B > 0 et xoi cs_B ^0 d'oii Xy cs_a ^ 9o cs_b < 1 - Xh cs_a 

e e e e e 

yci csjB > 0 et xgics_b ^ 0 d'oii yv cs^ ^ Qo cs-b ^ 1 - ye cs-a 

H convient de noter que les coefficients d'equilibre du constituant (I) ainsi etablis ne sont 
pas formellement equivalents a ceux que Ton obtiendrait en appliquant les jeux 
d'equations (6) h. (11) k la chaine de separation CS_B, moyennant la connaissance des 

e e e e 

masses molaires {MMgg, MMoo» MMog> MMgo}cs_b» ces masses molaires devant 
satisfaire, vis a vis des masses molaires des constituants (V), (I), (H), les inegalites decrites 
dans le paragraphe d). 

On peut, si besoin est, reviser le choix des masses molaires des constituants (V), (I), (H), 
en tenant compte des masses molaires des phases issues de la separation CS__B dans 
I'ecriture des inegalites du paragraphe d). 

j) A titre d^exemple, la demarche non limitative de realisation suivie pour obtenir la 
representation BRO pour le cas presente en illustration a consiste en differentes 

etapes: 



1 er depot 

,23 

generer avec un simulateur thermodynamique tel que par exemple le simulateur dit 
TVT PACKAGE') un ensemble d'etats 'e^' par une operation de depressurisation k 
volume constant a la temperature du reservoir (93.3 "^C), a partir de la pression de 
saturation du fluide (237.4 bar) jusqu'a 69 bar, en utilisant une equation d'etat de 
Peng-Robinson, d'emploi tres largement repandu dans le contexte petroiier, pour le 
fluide detaille k 16 constituants dont la composition initiale et les parametres sont 
donnes dans le Tableau 1 ci-apres : 



Tableau 1: Cas SPE3 - Description de reference du fluide de reservoir 



Nom 


Fraction 
molaire 


MM 


Tc 

op 


Pc 

psi 


Vc 

It /lom 


(0 




K (*) 
Okk' 


N2 


CI 


C02 


C2 


C3 


N2 


.0194 


28.016 


-232.33 


492.84 


1.4433 


0.0350 




N2 


0.000 


# 


# 


# 


# 


C1 


.6599 


16.043 


-115.76 


672-98 


1.5858 


0.01 30 




CI 


0.120 


0.000 


# 


# 




C02 


.0121 


44.010 


88.07 


1069.7 


1.5057 


0.2250 




C02 


0.000 


0.15000 


0.000 


# 


# 


02 


.0869 


30.070 


90.34 


709.67 


2.3707 


0.1050 




C2 


0.120 


0.00000 


0.150 


aooo 


# 


03 


.0591 


44.097 


206.28 


617.28 


3.2037 


0.1520 




C3 


0.120 


0.00000 


0.150 


0.000 


0.000 


104 


.0239 


58.124 


274.96 


528.95 


4.2129 


0.1918 




IC4 


0.120 


0.00000 


0.150 


0.000 


0-000 


04 


.0278 


58.124 


305.64 


550.56 


4.1007 


0.2010 




C4 


0.120 


0.00000 


0.150 


0.000 


0.000 


ICS 


.0167 


72.151 


370.11 


483.41 


4.9337 


0.2060 




ICS 


0.120 


0.00000 


0.150 


0.000 


0.000 


G5 


.0112 


72.151 


385.92 


489.36 


4.9817 


0.2520 




C5 


0.120 


o.ooo'oo 


0.150 


0.000 


0.000 


06 


.0181 


84.000 


452.25 


468.47 


5.9268 


0.2809 




C6 


0.120 


0.03000 


0,150 


0.010 


0.010 


C7 


.0144 


100.205 


512.64 


396.68 


6.9200 


0.3520 




C7 . 


.0,120. 


0.03000 


.0.150 


0.010 


0.010 


08 


.0150 


114.232 


564.64 


362.02 


7.881 1 


0.3992 




C8 


0.120 


0.03200 


0.150 


0.010 


0.010 


09 


.0105 


128.259 


613.31 


331,99 


8.7621 


0.4439 




C9 


0.120 


0.03400 


0.150 


0.010 


0.010 


CIO 


.0073 


142.286 


655.07 


304.15 


9.6591 


0.4869 




CIO 


0.120 


0.03600 


0.150 


0.010 


0.010 


C11. 


.0049 


147.000 


699.91 


349.54 


10.172 


0.4770 




C11 


0.120 


0.04300 


0.150 


0.010 


0.010 


C12P 


.0138 


210.000 


797.52 


324.79 


14.000 


0.5266 




C12P 


0.120 


0.04408 


0.150 


0.010 


0.010 
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Dans ee tableau, les parametres d'interactions binaires forment une matrice symetrique de 
trace nulle ; les cases ou figurent un symbole # sont remplies par symetrie. Les colotines 
manquantes doivent etre remplies d'une part par la propriete de symetiie, d' autre part par 
des parametres d'interactions binaires nuls. 

5 - A divers paliers de depressurisation repartis entre la pression de saturation du fluide 
initial et 69 bar, on a effectue une operation de vaporisation a pression constante (la 
pression du palier considere) de la phase huile de Tetat 'ea* considere par un gaz dont 
la composition est tenue pour representative du gaz d'injection en cours de simulation 
(dans le cas SPE3, la composition du gaz injecte varie en cours de simulation puisque 
10 c*est le gaz produit qui est reinjecte), chaque operation de vaporisation comportant 

plusieurs paliers de quantites de gaz injecte pour aller jusqu'a la vaporisation coinplete 
de la phase liquide. A chaque palier de quantite de gaz injecte ou les phases sont 
saturees correspond un etat 'edv'. 

Toujours avec le meme simulateur, les operations de ^6paration des phases huile et gaz 
15 de chacun des etats *ed' et 'Cdv' formant 1' ensemble « E » ont ete simulees, la chaine de 

separation CS_A utilisee comportant un premier etage a 26.7 ""C et 56.2 bar, un second 
etage de detente de la phase liquide issue du premier etage de separation a 26.7 °C et 
4.5 bar, une demiere detente de la phase liquide issue du deuxidme etage de separation 
a 15.6 ""C et 1,01325 bar, les phases gazeuses issues des differentes separations etant 
20 melangees. Comme le cas SPE3 presente une modification, au cours de la simulation, 

des conditions de pression du premier etage de separation (la pression passant de 56.2 
bar a 21.7 bar), les operations de separation des phases huile et gaz de chacun des etats 
'e' ont egalement ete simulees avec cette autre chaine de separation CS„B- 

II convient de noter que renchamement des di verses operations de simulation ci-dessus est 
25 aisement realisable avec la plupart des logiciels de simulations thermodynamiques 
disponibles pour Tindustrie. 

Les conditions necessaires decrites dans le paragraphe d) trouvees a partir des masses 
mol aires des produits de separation des phases huile et gaz de chacun des etats de 
r ensemble « E » par les chaines CS_A et CS_B sont en grammes/mole: 
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MMv< 21.592 



27.865 < MMi < 72,334 



MMh > 162.06 
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Les masses molaires attribuees aux trois constituants dans la representation BRO sent 
donnees dans le Tableau 2. 

1-es compositions a trois constituants des phases de I'ensemble « E » par les calculs decrits 
dans les paragraphes f) et g) ont ete obtenues en considerant les produits de separation de 
la chaine CS_A. La composition du fluide initial dans la representation BRO est donnee 
dans le Tableau 2 ci-apres : 



Tableau 2: Gas SPE3 - Description BRO du fiuide de reservoir 


. Nom du 
constituant 


MM 


Fraction 
molaire 


V 


21.592 


.856714 


I 


62.334 


-092358 


H 


162.06 


.050928 



Proprietes des phases , * 

La regie de Gibbs, qui donne le nombre de degres de liberte d'un systeme 
thermodynamique a partir du nombre de constituants et du nombre de phases, appli^uee a 
un melange a trois constituants, donne trois degres de liberte lorsque deux phases 
hydrocarbure sont a I'equilibre : les proprietes de phases dependent de la pression, la 
temperature et d'une variable compositionnelle. 

Pour chaque phase fluide pour laquelle on connait, pour des conditions de pression et 
temperature donnees, les valeurs de proprietes de phase utiles pour la modelisation des 
ecoulements (facteur de compressibilite, viscosite, ..), on pent rentrer ces proprietes de 
phase sous forme de tables fonctions de la pression, de la temperature/ et d'un index de la 
composition dans la representation a trois constituants. Ces donnfes peuvent etre utilisees 
ensuite comme donnees d'entrees jdans un simulateur, le simulateur disposant en general 
de methodes internes d'interpolation ou d'extrapolation pour estimer les proprietes de 
phases en des points intermediaires. / 

k) Altemativement a T usage de tables, on peut rechercher, par exemple par 
regression, les parametres de coiTelations ou d'equations d'etat permettant de calculer, par 
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zone thermodynamique, les proprietes des phases en fonction de la composition et des 
conditions de pression et de temperature, ces parametres etant ensuite introduits dans le 
modele de simulation en lieu et place des tables precedemment citees. 

Une situation frequente est de disposer d'une equation d'etat fiable pour le fluide de 
5 reservoir decrit de fagon detaillee dans la representation de base. 

Ainsi, lorsqu'une equation d'etat (et une correlation associee pour le calcul des viscosites) 
definie pour la representation detaillee \ Nrb constituants, permet de niodeliser le 
comportement thermodynamique des fluides au cours des etapes de pressions et 
temperatures d'un chemin thermodynamique suivi dans une ou plusieurs zones 

10 thermodynamiques, on peut utiliser cette equation d'etat pour generer 1' ensemble de 
reference de compositions de fluides et les proprietes deja citees des phases 
correspondantes, en y adjoignant les parametres des phases dans I'equation d'etat, tels que 
le CO- volume, le terme d' attraction, la pression critique, qui peuvent etre consideres 
conrmie des proprietes specifiques des phases. Pour chacune des compositions de fluides 

15 obtenues, cette meme equation d'etat, ou une autre, peut etre utilisee pour generer les 
compositions et proprietes des melanges fluides produits par une chaine de separation, et 
permettre ainsi 1' application de la methode de pseudoisation decrite dans les paragraphes 
a) a i). 

Lorsque les parametres des pseudo-composants dans une equation d'etat, et dans une 
20 correlation pour le calcul de viscosite, sont ajustes, par exemple par regression, pour 
reproduire les proprietes specifiques des phases telles que le co-volume, le terme 
d' attraction, la temperature critique, .. obtenues avec la representation detaillee, on obtient 
une equation d'etat et une correlation utilisable pour les calculs de proprietes de phase 
telles que le facteur de compressibilite, la viscosite, mais cette premiere Equation d'etat, 
25 designee ci-apres en abrege par EOS_PRO, n'est pas forcement valide pour les calculs 
d'equilibre notamment si le respect d'une condition d'egalite des fugacites a I'equilibre n'a 
pas ete pris en compte. 

Une demarche non limitative de realisation pour obtenir les parametres de I'equation d'etat 
EOS__PRO de la representation BRO est celle suivie pour le cas 'SPE3' presente en 
30 illustration. L'equation d'etat utilisee, celle de Peng-Robinson, est rappelee ci-apres pour 
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aider h, la comprehension de la demarche, mais la demarche est aisement transposable ^ 
d'autres equations d'etat. 

1) L'equation de PengrRobinson s'ecrit pour la phase P : 

(15) p = ^- — 5^9^ 

Vp-bp (Vp + bp) -2bp 

ou R est la constante des gaz parfaits, T la temperature, P la pression, V le volume molaire 
de la phase P, b son co-volume et a(T) son terme d' attraction, ces demidres quantites etant 
exprimees par : 

(16) bprrZcibi 



(17) a(T)p = 2:cj 



5:ciaij(T) 



avec aij(T) = ^J^) (1 - 5ij) 



10 ou Ci designe la fraction molaire du constituant i dans la phase P, fraction inolaire dont les 
d6pendances physiques sont omises pour alleger la presentation, b\ et ai(T) respectivement 
le co-volume et le terme d'attraction du constituant i, 5ij le parametre d'interaction binaire 
entre les constituants i et j, les formant une matrice symetrique de trace nulle. 

Sous forme adimensionnelle, 1 'equation de Peng-Robinson s'exprime comme: 

? 

3 2 2 2 3 

15 (18) Zp -(l-Bp)Zp +(Ap-3Bp - 2 Bp) Zp- (ApBp-Bp -Bp) = 0 

ou Z est le facteur de compressibilite de la phase P, et B et A sont les formes 
adimensionnelles des co-volume et terme d' attraction de la phase P (les dependances des 
variables avec la pressiOn et temperature etant orhises pour alleger Tecriture des 
equations): 

20 (19) Bp = bp^ 



(20) Ap = ap —5—7 
R T 
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m) Le facteur de compressibilite de la phase P est uniquement fonction des deux 
parametres de phase adimensionnels A et B. II s'ensuit qu'estimer les parametres par 
constituant dans la representation BRO (a savoir bv, bi, bH, av(T), ai(T), aH(T), 5vi, 6vh, 
8ih) permettant de reproduire les valeurs de a(T) et b, obtenues pour chaque phase aux 
differents paliers des simulations PVT avec la representation de base (detaillee h. 16 
constituants), permet de reproduire les valeurs de reference des facteurs de compressibilite 
des phases. La procedure de regression utilisee pour obtenir les differents parametres de la 
representation BRO est la minimisation de fonctions objectif de la forme: 



(21) 0(p) = Z [<Pg(p) - og] + Z ^o [(Po(p) - oo] 



10 ou la quantite observable o par phase et par etat e, ici obtenue avec la representation 
detaillee a 16 constituants, est modelisee par une fonctionnelle (p parametree par les 

composantes pn du vecteur p, et ou les quantites tu sont des poids attribues aux differentes 
observables. 

Le minimum de 0(p) est obtenu par la resolution du systeme d'equations, de meme 

15 dimension que le vecteur p, forme par la nullite des derivees partielles de la fonction 0(p) 
par rapport aux differents parametres p^ : 

(22) ap^O(p) = 0 



soit: 



(23) S [<Pg(p) - og] 3p^(pG(p) + S ^o[(poG) - Oo] 5p^<Po(p) = 0 

20 n) Pour la recherche des valeurs des parametres bx (K=V,I,H) de la representation 
BRO, la fonction objectif utilisee dans le cas SPE3 est done: 



(24) 0(bv,bi,bH) = S ^g 



SyK^bK- be 

K 



+ S ^O 



E^K^bK- bo 

K 
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Le sy Sterne a resoudre est simplement un systeme lineaire de trois equations k trois 
inconnues: 



(25) 



nivv bv + nivi bi + nivH — 
mp/ by + mu bi + mm hn = rj 
niHv bv + niHi bi + niHH bn = th 



les termes hikk- tk des seconds membres des .equations etant donnes par : 

eeE ^ ^ e€E . ^ ^ ^ e 

5 niKK- = Z (^G yK^ Yk'^ + y^K ^k) ; Tk = S ( bo + vdg yK be) 



o) Pour rechercher les parametres av(T), ai(T), aH(T), 5vi, 5vh, 5ih par simplement des 
resolutions de systemes lineaires de trois equations a trois inconnues, Tequation (17) peut 
etre reecrite comme la somme de deux termes, le second terme comprenant les interactions 
binaires : 

10 (26) a(T)p = a,(T)p%a2(T)p 
avec : 

ill) a,(T)p = EciVii(T) et 

i ■ 

a2(T)p = --£ ScjCi aijou cXij(T) = Si| 

j i 

Dans le terme a2(T)p, les parametres 0Cij(T) forment un ensemble de trois parametres 
15 independants avi(T), avH(T), aiH(T). 

En utilisant les tefmes ai(T), obtenus lors des simulations PVT avec la representation 
detaillee a 16 constituants, comme quantites obser\'ables o par phase et par etat e, les 



inconnues etant les quantites '\/av(T), -^/aKT), VaH(T), le probleme a resoudre est bien 
analogue & celui decrit pour I'obtention des parametres de co-volumes des constituants de 
20 la representation BRO. De meme, les parametres oCij(T) sont obtenus en resqlvant le 
systeme de trois equations a trois inconnues forme & partir des observables a2(T). 
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p) Les parametres d' interactions binaires peuvent etre ensuite obtenus par: 
(28) 

5 ^ avi(T) _ avH(T) amCT) 

q) Certains simulateurs acceptent, comme parametres des constituants pour T equation 
d'etat, rintroduction directe des co-volumes bj des constituants et de la matrice des 
parametres ay(T) intervenant dans le calcul du terme d'attraction (cf. I'equation (17)), cette 
matrice devant 8tre definie le cas 6cheant pour differentes temperatures. 

Au cas ou le cas traite necessiterait une dependance en temperature (temperature variable 
dans le reservoir), la demarche prec^dente pourrait etre adaptee, par exemple en generant 
Tensemble «E» par des operations PVT k diverses temperatures de reservoir, en 
considerant des etats e a diverses temperatures pour obtenir les co- volumes par constituant 
(ceux-ci devant etre les memes pour les differentes temperatures), en considerant les etats 
e obtenus a une temperature donnee pour generer la matrice des parametres aij(T^) 
correspondante et en repetant la demarche pour autant de temperatures que necessaire. 

r) Lorsque le simulateur n'accepte pas Tintroduction directe du co-volurne par 
constituant hi et des matrices de parametres aij(T^): 

- le co-volume par constituant est en general obtenu a partir de rintroduction des 
parametres Tci, Pci, optionnellement Qbi (par defaut tous egaux a 0.077796) par : 

(29) bi=:RQbi|| 

- les quantites ai(T) sont en general obtenues par: 

(30) ^RT) = I 1 + Ci f 1 - ^ /^ll avec Ci = S C„ co^" 



n=0. 



La quantite aOi est en general obtenue a partir de rintroduction des parametres Tcj, Pci, 
optionnellement £2ai (par defaut tous egaux a 0.457235) par: 
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(31) aOi = R^aai^ 

Dans le polyn6me Q de degr6 dp, les parametres oh sont les facteurs acentriques des 
constituants en general consideres comme constants, le degre du polyn6me et les 
constantes Cn pouvant varier d'un simulateur a T autre. 

5 s) Lorsque les pressions, temperatures et volumes critiques par constituant sont 
requis, on pent utiliser^conmie observables les pressions, temperatures et volumes critiques 
des phases de Tensemble «E», definies par des lois de melanges identiques h celles 
utilisees pour les co-volumes : 

(32) Pcp = XciPci ; Tcp=:XciTci ; Vcp = SciVci 

i i i 

10 L'obtention respectivement des parametres Pcv, Pci, Pch, Tcv, Tci, Tch et Vcy, Vci, Vch 
peut done s'effectuer par resolution successive de trois systemes de trois equations a trois 
inconnues totalement analogues a ceux generes pour le calcul des, co- volumes par 
constituant, developpe precedemment. 

t) Les parametres Qb^ (K=V,I,H) sont alors obtenus, a partir de 1' equation (29) et des 
15 parametres bK, Pck et Tck tels que detemiines dans les paragraphes n) et s), par: 

(33) ^2bK = f^ ■ 

u) Dans Tequation (30), dans le cas le plus frequent ou la temperature de reservoir est 
consideree comme constante, on peut s'affranchir de la dependance en temperature des 



quantites ^/av(T), V^iCT), V^h(T), en adoptant une valeur nulle pour les facteurs 
20 acentriques des constituants V, I, H, d'ou: 

(34) VaK(T) ^ , VK, K=V,I,H 

et on obtient alors les paiameties iiaK (K=VJ,H) , a partir des parametres ^Ja^iT), Pck et 
Tck tels que determines dans les paragraphes o) et s), par: 

/'JC^ r-i Pck aK(T) 

(35) OaK= R^TCK^ 
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Une autre possibilite, utilisee dans le cas SPE3, est T attribution aux parametres i2aK 
(K=V,I,H) de la valeur par defaut, Dans ce cas, les parametres ^JaQ^ sont determines avec 
les Pck et Tcr tels que determines dans le paragraphe s), et les parametres sont ensuite 
obtenus par rearrangement de F equation (30): 



(3«) Ck = ^ 



^ \frcK 

puis, les facteurs acentriques cok ont ete recherches conrnie solutions des trois equations 
polynomiales (K=V,I,H): 

(37) ^Cn(OK -Ck = 0 

n=0 

v) Lorsque pour le cas considere, la temperature du reservoir est variable, on peut 
10 reecrire I'equation (30) comme : 

(38) VaK(n-PK + YK'nK(T^) 

avec les quantites : tIk(T') = 1 - ^ calculees en utilisant les temperatures critiques 

des constituants telles que determinees dans le paragraphe s). 

On peut ensuite rechercher les quantites Pk et Yk Ok = V^Ok et Yk = V^Ok Ck) comme des 
15 parametres independants par regression lineaire des couples de valeurs tikCI*), '\/aK(T^), 
ces demieres valeurs etant telles que determinees dans les paragraphes o) et q). 

Les parametres QaK et Ck sont alors obtenus par: 

(39) aaK = p^ et = g 

Les facteurs acentriques recherches sont solutions des trois equations polynomiales 
20 (K=V,I,H): 
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(40) SCnO)K"-l; = 0 

w) Des facteurs de correction de volume sont souvent utilis6s pour ajuster les densites 
des phases. Lorsque des facteurs de correction de volume par constituant sont requis, on 
peut utiliser comme observables les corrections de volume des phases de Tensemble « E », 
5 definies par des lois de melanges identiques a celles utilis^es pour les co- volumes: 

(41) Cvp=2ciCvi 

i 

L'obtention des parametres Cvv, Cvi, Cvh, peut s'effectuer alors par resolution d'un 
systeme de trois equations a trois inconnues totalement analogue a celui genere pour le 
calcul des co-volumes par constituant developpe dans le paiagraphe n). 

10 x) Dans le cas SPE3, les calculs de viscosite des phases hydrocarbure ont ete effectues 
avec une correlation largement utilisee dans le domaine de Tingenierie p6troliere, pelle de 
Lohrenz, John, Bray, B.G. et Clark, C.R.: « Calculating Viscosities of Reservoir Fluids 
from Their Compositions » , Journal of Petroleum Technology, 1964, pp. 1171-1176. Les 
parametres de phases utilises explicitement dans cette correlation sont les pressions, 

15 temp6ratures, volumes critiques des phases. La determination des parametres critiques 
(pression, temperature, volume) des constituants V, I, H, telle qu'exposee dans le 
paragraphe s), a permis une modelisation appropriee des viscosites des phases avec la 
representation BRO. La demarche est aisement transposable a d'autres correlations, 

y) La ligne directrice suivie dans Tensemble des paragraphes n) a w).est de rechercher 
. 20 les parametres d'equation d'etat EOS_PRO utiles aux calculs, dans les conditions de 
reservoir, de densites des phases (via les facteurs de compressibilite), et de viscosites par 
des resolutions simples (non iteratives) de systemes d'equations lineaires. 

Son efficacite pour reproduire les co-volumes et termes d' attraction des phases est 
illustree, pour les fluides du cas SPE3, par les Fig. 6 et 7^ pour une operation de 
25 depressurisation a volume constant a la temperature de reservoir' (93.33 a partir de la 
pression de rosee bars du fluide de reservoir initial (237,4 bars) jusqu'a 69 bars. Sur 
chacune des figures, sont illustrees les quantites calculees a partir de Fequation d'etat de la 
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representation detaillee de reference et les quantites calculees a partir de Tequation d'etat 
EOS_PRO de la modelisation BRO, ces quantites etant pratiquement confondues. 

Cette meme ligne a ete utilisee avec le meme succes pour traiter d'autres cas non presentes 
ici. On a pu constater qu'une telle ligne est largement adaptable a des contextes de 
5 temperature variable. Par ailleurs, les traitements necessaires ne representent aucune 
difficulte particuliere de codage et ont ete aisement automatisables. 

Les parametres de rEOS_PRO utilises pour le cas SPE3 en illustration sont donnes dans le 
tableau 3 ci-dessous : 



Tableau 3: Cas SPE3 - Parametres de i'EOS^PRO de la representation BRO 


Nom 


Tc 


Pc 

psi 


ftMbm 


CO 




Cv 
ftMbm 






V 


I 


H 


V 


-61.212 


666,02 


1.869 


-0.02043 


0.079345 


-0.0686 




V 


0.0 


0.058447 


0.002065 


1 


368.292 


581.21 


4.389 


0.31255 


0.082877 


-0.0621 




1 


# 


0.000000 


0.001604 


H 


663.138 


263.46 


10.99 


0.21233 


0.065003 


0.2664 




H 


# 


# 


0.000000 



10 

z) Dans le cas SPE3, on a choisi d'utiliser pour I'EOS^PRO une equation d'etat du 
meme type que celle utilisee pour la representation de base. Ceci n'est pas obligatoire. II 
est possible d'utiliser une equation d'etat d'un autre type, des lors que cette autre equation 
d'etat utilise des parametres de phase identifiables aux parametres de phase utilisees par 
15 Tequation d'etat de reference. 

Plus generalement encore, on pent obtenir les parametres necessaires a la representation 
BRO par diverses techniques dites de regression, inversion, optimisation, que Ton peut 
trouver dans une litterature abondante. On peut par exemple les rechercher par une ou 
plusieurs resolutions iteratives de systemes d' equations formes par des derivees partielles, 
20 par rapport a des sous-ensembles ou a I'ensemble des parametres a la representation BRO, 
d'une fonction objectif du type de Tequation (21) ou d'un autre type. 

On pourrait meme chercher a utiliser un formalisme de ce type pour optimiser Tune ou 
plusieurs des masses molaires des constituants V, I, H, en la (les) faisant apparaitre 
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explicitement, - par remplacement dans la fonction objectif des compositions des phases 
par leurs expressions donnees par les equations (8) et (9) et en choisissant des 
observables sensibles a leur(s) valeur(s). Un telle demarche pourrait etre envisageable si 
Ton se trouvait dans une situation, non rencontree jusqu'a present, ou la ligne directrice 
5 suivie dans I'ensemble des paragraphes n) a w) conduisait a une determination de 
parametres de la representation BRO ne permettant pas une reproduction fidele des 
observables. 

Enfin, pour des raisons connues en soi, on peut etre amene a donner a priori des valeurs & 
certains des parametres de Tequation EOS_PRO et n'effectuer alors des calculs 
10 d'estimation que pour un nombre reduit de parametres. 

aa) La premiere expression generalement attendue des resultats de simulation sont des 
resultats en volume. Les observables que Ton doit alors chercher a restituer sont les 
volumes molaires des phases en sortie de s6paration. Pour le cas SPE3, on a utilise deux- 
correlations distinctes. Tune pour le calcul des volumes molaires des phases gaz en^^ortie 
15 de separation, T autre pour le calcul molaire des phases hiiiles en sortie de separation, les 
conditions de sortie de separation etant les conditions standard usuelles de la profession & 
savoir une atmosphere et 15.56 ""C (60 °F). 

De nombreuses corr61ations ou equations d'etat sont disponibles pour estimer le volume 
molaire d*une phase gaz notamment dans des conditions proches des conditions 

20 atmosph6riques. A litre indicatif, c'est la correlation de P.M. Dranchuk et LH. Abou- 
Kassem (presentee dans le papier : « Calculating Z Factors For Natural Gases Using 
Equations of State », JCPT, Juillet-Septembre 1975) qui a ete utilisee pour le cas SPE3 oii 
les parametres utiles au calcul du volume molaire du gaz sont les pressions et temperatures 
critiques des constituants (V) et (I), parametres qui ont ete regresses de fa9on iterative de 

25 fagon a reproduire les volumes molaires des phases gaz obtenus par les simulations PVT 
de separation avec la representation a 16 constituants. 

La variation, en fonction de la fraction molaire du constituant I, du volume molaire des 
phases huile en des conditions proches des conditions atmbspheriques a ete trouvee, dans 
le cas SPE3, pratiquement lineaire, plus generalement, pour des cas non presentes ici, 
30 representable par un polynome de degre peu eleve (degre deux a trois). Pour le cas SPE3, 
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c'est une dependance lineaire qui a ete adoptee dont les parametres ont ete obtenus par 
simple regression lineaire. 

Les parametres necessaires aux calculs des volumes molaires en conditions de separation 
ont ete obtenus pour les deux chaines de separation CS_A et CS_B avec la meme methode, 

Equilibres thermodynamiques 

bb) Les calculs d'equilibre peuvent etre effectues soit a partir de donnees de 
coefficients d'equilibre, soit ^ partir de Tecriture de I'egalite des fugacites de tous les 
constituants dans chaque phase, ces deux possibilites etant examinees respectivement ce 
qui suit. 

cc) Selon la regie de Gibbs, pour le calcul des equilibres dans le reservoir, les 
coefficients d'equilibre' des pseudo-composants (V), (I), (H) tels que calcules par le jeu 
d' equations (10) peuvent etre introduits dans des tables comme fonctions de la pression, de 
la temperature, et d'un index de la composition dans la pseudo representation, Ces donnees 
peuvent etre utilisees ensuite comme donnees d' entrees dans un simulateur, le simulateur 
disposant en general de methodes internes dMnterpolation ou d' extrapolation pour estimer 
les coefficients d'equilibre en des points intermedi aires. On pent aussi, avant introduction 
des donnees dans un simulateur, ajouter des points aux tables en utilisant des methodes, 
connues en soi, d' interpolation ou d' extrapolation ou utiliser, a la place des tables, des 
correlations baties pour restituer les coefficients d'equilibre. 

Une dependance des coefficients d'equilibre avec, tout a la fois la pression, la temperature 
et rindex de composition, n'est pas toujours necessaire. Ainsi, dans nombre de cas 
d'application, la temperature du reservoir est consideree comme constante. Dans certains 
cas, on peut etre amene a abandonner la dependance par rapport a un index 
compositionneL Ces cas ou Ton reduirait le nombre des degres de liberte des coefficients 
d'equilibre restent dans le domaine de f utilisation de la methode proposee. 

dd) Pour le calcul des equilibres a Tissue d'une chame de separation, les coefficients 
d'equilibre des pseudo-composants (V) et (H) sont respectivement infini et nul quelques 
soient les conditions de pression et de temperature en sortie de chaine de separation et 
quelque soit la composition du melange soumis a la separation. Si besoin est, dans un 
simulateur qui n'accepterait pas le concept de coefficients d'equilibre infini ou nul (ce 
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concept etant sou vent implemente en forgant T absence du constituant respectivement dans 
la phase huile et dans la phase gaz), on pent utiliser une valeur finie tres grande et une 
valeur finie tres faible, toutes deux positives, en lieu et place de coefficients d'equilibre 
infini ou nuL Pour les donnees relatives aux coefficients d'equilibre du constituant (I) 
5 necessaires a un calcul d'equilibre, le jeu d'equations (11) indique une definition diff6rente 
des coefficients d'equilibre du constituant (I) selon que Ton considere la separation d*une 
phase 'gaz' issue du reservoir, ou la separation d'une phase *huile' issue du reservoir. Les 
entrees, en principe necessaires aux calculs d'equilibre, sont done, par phase, une table de 
coefficients d'equilibre du constituant (I) en fonction d'un index compositionnel de la 
10 phase soumise a la separation, de la pression et de la temperature en sortie de la chmne de 
separation considered Les plages de variation des compositions des phases huile et gaz 
etant en general distinctes, il peut etre pratique de 'fusionner' en une seule table les 
donnees utiles pour les phases huile et gaz, notamment si le simulateur n'accepte pas 
r introduction de donnees distinctes pour chaque phase, 

15 Ces donnees sont utilisables cpmme donnees d' entrees dans un simulateur, le' simulateur 
disposant en general de methodes internes d' interpolation ou d' extrapolation pour estimer 
les coefficients d'equilibre en des points intermedi aires. On peut aussi, avant introduction 
des donnees dans un simulateur, ajouter des points aux tables en utilisant des methodes, 
connues en soi, d'interpolation ou d' extrapolation ou utiliser, a la place des tables, des 

20 correlations baties pour restituer les coefficients d'equilibre. 

ee) Si I'equation d'etat EOS_PRO est batie pour restituer les proprietes des phases, elle 
n'est pas forcement valide pour les calculs d'equilibre, c'est-a-dire pour permettre de 
restituer les compositions des phases de I'ensemble «E » obtenues dans le paragraphe g). 
Comme precedemment, on utilise ci-apres I'equation d'etat de Peng-Robinson comme 
25 support de presentation. 

Pour cette equation d'etat, les coefficients d'equilibre, d'apres la condition d'egalite des 
fugacites de chaque constituant dans chaque phase a I'equilibre, s'obtiennent par: 
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N 

(42) Ln(Ki) = Do + Di bi + ^^2^ X (fe ~ ''V ^'^^^ 



avec les parametres adimensionnels de phase definis selon les equations (19) et (20) et: 



(43) 



Zp + (l-V2)Bp^ 



^Zo '-' Bo J 
P "Zo-1 Zg - 1 AqLq 



Lp = Ln 
< Do = Ln| 

°^ = Rt[ Bo " Bg "^2V2Bo''2^/2Bg^ 

et les facteurs de compressibilite des phases calcules a partir de T equation (18). 

Pour traiter un cas ou la temperature de reservoir est constante, le probleme est de 
determiner les neuf parametres independants formant un vecteur p : by, bi, bn, aw, an, ann, 
avi, avH, aiH, par example par minimisation d'une fonction objectif construite en utilisant 
comme observables ok les coefficients d'equilibre de I'ensemble «E» tels que calcules 
par les equations (10). La fonction objectif selon Tequation (21) peut alors etre adaptee 
comme: 



eeE 



10 (44) 0(p) =5: S [ (Pk(p) - Ok] 



e K=V.I,H 



15 



OU les fonctions (Pk(p) sont les exponentielles des seconds membres des equations (42), 

e 

Altemativement, dans T equation (44), on pourrait utiliser comme observables Ok les 

e — > 

logarithmes des coefficients d'equilibre de Tensemble «E» et comme fonctions (Pk(p) 
sont les seconds membres des equations (42). 

En suivant la demarche de K.E. Starling dans le papier: «A New Approach for 
Determining Equation-of-State Parameters Using Phase Equilibria Data » , SPE 1481, 

e 

SPEJ, Decembre 1966, on peut choisir les poids Wk comme les inverses des valeurs des 

e 

observables Ok- 



1 er depot 



39 . 

La nullite des derivees partielles de la fonction O(^) par rapport aux neuf parametres 
d'interet foumit un systeme de neuf equations non lineaires a neuf inconnues, resoluble 
uniquement de fa^on iterative, de nombreuses methodes adaptees a la resolution des 
systemes non lineaires etant discutees dans la litterature. 

Un calcul iteratif de ce type necessite une premiere approximation du vecteur de 
parametres. On a pu constater, notamment pour le cas SPE3, qu'il est opportun de titer 
profit des parametres obtenus pour rEOS_PRO^pour demarrer les calculs iteratifs. 

Plus generalement, notamment pour reproduire des coefficients d'equilibre a diverses 
temperature, on peut considerer un autre ensemble ou sous-ensemble de parametres pour le 

vecteur p et d'autres techniques de recherche de parametres. 

ff) A titre d' illustration, pour le cas SPE3, ce sont les six parametres Pck et Qb^ qui 

— > 

ont constitu6 le vecteur p, les valeurs des autres parametres etant reprises de Tequation ^ 
EOS_PRO. La procedure utilisee dans le cas SPE3 est une procedure connue en ^soi en 
cours d'amelioration. L'ensemble des parametres de I'EOS^EQ utilises pour le cas|;SPE3 
est donne dans le tableau 4 suivant : , 



Tableau 4: Cas SPE3 - Parametres de I'EpS^EQ de !a representation BRO 


Nom 


op 


psi 


CD 








V 


I 


H 


V 


-61.212 


848.47 


-0.02043 


0.089345 




V 


0.0 


0-058447 


0.002065 


1 


368.292 


274.23 


0.31255 


0.070000 




I 


# 


o.oooodo 


0.001604 


H 


663.138 


309.90 


0.21233 


0.065003 




H 


# 




O.OOOOOO 



gg) Dans le cas SPE3, les donnees introduites pour les calculs d'equilibre en surface 
pour chacune des chaines de separation CS_ A et CSJB sont les coefficients d'equilibre du 
constituant (I) sous forme d'une table unique pour la separation des phases huile et gaz, 
table obtenue comme decrit dans le paragraphe dd) avec pour seiil abscisse Tindex de la 
composition de la phase avant separation. 
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hh) Les Fig. 2 a 5 comparent les resultats, sur les quinze annees d'exploitation simulees 
du cas SPE3, d'une simulation compositionnelle a trois pseudo-composants BRO (resultats 
designes par 'temaire'), ou la representation compositionnelle a ete obtenue par 
application de la methode selon I'invention, a ceux d'une simulation compositionnelle 
5 detaillee a seize constituants de base (resultats designes par 'reference'). 

La simulation compositionnelle de reference utilise une meme equation d'etat de Peng- 
Robinson pour les calculs d'equilibre et pour les calculs de proprietes de phase et pour 
toutes les zones thermodynamiques. La simulation BRO a trois constituants utilise, pour 
les calculs dans les mailles de reservoir, I'equation d'etat (EOS_EQ) pour les calculs 
10 d'equilibre et Tequation d'etat (EOS^PRO) pour les calculs de proprietes de phases. Pour 
les calculs en conditions de surface, la simulation BRO utilise des coefficients d'equilibre 
introduits sous forme de tables et des correlations simples pour le calcul des volumes des 
phases. 

Les resultats presentes Fig. 2, qui representent revolution de pression en fond de puits 
15 producteur, et Fig. 5, qui representent revolution de saturation en huile (ou condensat) 
dans la maille perforce la plus profonde a I'endroit du producteur, sont panni les resultats 
les plus revelateurs du bon calcul des equilibres et des proprietes de phases dans les 
conditions in situ. 

On peut constater que la modelisation BRO foumit une solution en pression pratiquement 
20 confondue avec la solution de reference (a 16 constituants) et une solution d' evolution de 
saturation en huile en bon accord avec la solution de reference sur tout le temps simule. La 
maille utilisee en Fig. 5 est la meme qiie celle utilisee dans le papier de reference de D.E. 
Kenyon et G. Alda Behie et dans le papier de W.H. Goldthorpe cite dans I'etat de la 
technique. On peut constater que les resultats de la Fig. 5, juges par comparaison avec les 
25 resultats obtenus par W,H. Goldthorpe avec une modelisation BO, montrent une bien 
meilleure qualite de reproduction de revolution de la saturation en huile dans le reservoir 
alors meme que la procedure utilisee pour I'obtention des parametres de TEOS^EQ est 
ameliorable. 

ii) Les resultats presentes Fig. 3 et 4 de debit d'huile et de production cumulee d'huile 
30 sur toute la duree du temps simule montre que la modelisation BRO foumit une solution 
tres proche de la solution de reference. 
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La reduction de temps de calcul par la modelisation BRO est considerable, le facteur 
d*acceleration etant de 46. 

Delumping 

jj) On repr6sente tout d'abord le reservoir sous la forme d'un reseau de mailles (m) 
5 dont chacune constitue un volume 616mentaire rempli de fluide(s) sous forme de une ou 
plusieurs phases, avec au moins une phase non aqueuse. Les phases non aqueuses restent 
dites phases hydrocarbure meme si elles peuvbnt contenir certains composants autres que 
des hydrocarbures tel que T azote, le dioxyde de carbone, Tanhydride sulfureux, 

kk) On definit, pour chaque zone ou domaine thermodynanndque, le fluide par une 
10 . representation d6taillee & Nrb composants et/ou pseudo-composants. n est a noter que Ton 
peut traiter les cas qui requierent I'usage de plusieurs representations thermodynamiques, 
par exemple si Ton peut distinguer plusieurs chemins thermodynamiques locaux durant la 
modelisation (on peut etre amene a distinguer des zones produites uniquement par 
depressurisation et des zones soumises a une injection de gaz). Plusieurs zones ou 
15 domaines de variations des grandeurs thermodynamiques ou de composition, designes 
souvent par les specialistes comme zones thermodynamiques, pourront done 8tre definis et 
utilises. 

11) Par zone thermodynamique pour laquelle on adopte une representation groupee des 
fluides, on determine une equation d'etat EOS_PRO construite, prealablement a la 
20 simulation dynamique de reservoir avec la representation groupee, pour reproduire les 
parametres de phases, dans Tequation d'etat de la representation detaillee, des fluides 
hydrocarbure au cours de chemins thermodynamiques tenus pour representatifs de ceux 
que vont suivre les fluides hydrocarbure pendant la simulation maillee, ' 

mm) On realise, de fa9on connue en soi, une simulation compositionnelle avec un 
25 nombre limite de constituants ou les proprietes des phases sont calculees par une equation 
d'etat EOS_PRO, ladite sirnulation permettant de calculer au moins dans chaque maille 
(m) et a des pas de temps consecutifs (t, t+1, etc.), une pression dans une phase 
hydrocarbure (pm), la temperature (Tm) (si celle-ci varie), les debits molaires par phase 
d'injection ou de production (SgJn) et (Sdn). pour chaque paire de mailles (m,h), les debits 
30 molaires des phases liquide (uomh ) et vapeur (Ug,lih), la fraction vapeur (0m). On memorise 
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ces quantites, ainsi que les paramdtres dans T equation d'etat EOS^PRO pour les phases 
hydrocarbure liquide et vapeur (typiquement les co-volumes et termes d' attraction) 
intervenant dans T expression des coefficients d'equilibre. Pour eviter de les recalculer, on 
peut memoriser egalement les facteurs de compressibilite des phases (Zom) et (7^^ 
5 obtenus par resolution de Tequation d'6tat EOS_PRO. 

Pour une equation d'etat de Peng-Robinson, en conservant les notations introduites par les 

equations (16), (17), (26) et (27), les parametres de phase sont designes par (bom) et (bcm) 

pour les co-volumes des phases, (aoL) et (acm), - ou (a,om), (aiom), (azom), (sligL) pour les 
termes d' attraction des phases, 

10 Les differents parametres utiles au calcul des coefficients d'equilibre, lorsque evoques 
d'une fagon g^nerale, seront designes ci-apres par (En^\ I'indice P etant pour chaque 
phase huile et gaz, et I'indice n un numero de reference pour un parametre particulier. 

nn) On estime au pas de temps t+1 la fraction molaire de chaque constituant i dans la 
composition detaillee globale (ztm) du fluide hydrocarbure dans la maille (m) en 
15 connaissant les nombres de moles de chaque constituant i dans la representation detaillee 
dans chaque phase au pas de temps t, respectivement (No^) pour la phase huile et (Ngi^) 
pour la phase gaz a partir des equations ci-apres: 



Nrb 



Nrb 



(45) 




N,l = Nom + Ng, 



>in 



(46) 



Xi, 




Ng4 




^irn N] 



m 



At (yim> Sgm + Xi, 



Som) - At X (Yim' Ugmh + ^im' Uomh) 
heJ(m) 



20 (47) 



t+1 



FT 



ni 



(48) 



n'™ = - At (SgL + Soln) - At X (Ugilh + u„U) 



heJ{m) 
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Les equations (45) et (46) permettent respectivement de calculer, au pas de temps (t) et 

t 

dans la maille (m), les nombres de moles par phase hydrocarbure et global Nm a partir de 
la connaissance des nombres de moles par constituant et par phase, et les fractions 
molaires des constituants par phase et globales< 

5 L' equation (48) traduit le bilan molaire total hydrocarbure sur la maille (m), en tenant 
compte des echanges de matiere, pendant Tintervalle de temps At, avec toutes les mailles 
(h) voisines de (m) qui forment I'ensemble J(n^). 

Dans Tequation (47), Tecriture des termes (y[|^) et (Xj^) , - dans lesquels m'=m pour un 
debit de la maille (m) vers la maille (h) ou dans le puits, et m'=h pour un debit de la maille 
10 (h) vers la maille (m), et m' cbrrespondant au fluide injecte dans le cas des puits 
d' injection, S 6tant alors negatif suppose implicitement rutilisation d'un schema simple 
amont pour les flux compositionnels. Une ecriture plus gen6rale de ces termes est 
(yimh) et (Ximh) OU ct yi,^ decrivent les compositions des phases liquide et gaz/ 
obtenues d'une fagon connue en soi, s'ecoulant entre les mailles (m) et (h). 

15 oo) On determine, au pas de temps (t+1) et dans chaque maille (m), les coefficients 
d'dquilibre (K^m) du constituant i dans la representation detaillee en utilisant une 
formulation ou apparaissent explicitement les parametres Ejcp, en y introduisant les 
parametres variables (Enp^) issus de la simulation avec la representation groupee, et, le 
cas echeant, en utilisant certaines des variables compositionnelles decrites dans le 

20 ' paragraphe precedent calculees au pas de temps (t). Ainsi, pour une equation d'etat de 
Peng-Robinson, les coefficients d'equilibre donnes par Tequation (42) peuveiit etre 
reecrits sous la forme: 

^ ■ _t+i 

t+i t+1 t+l i+l I t+l P m 

(49) Ln(Kim)=Dom +Dim bj + D2m 'Vai(T„,) + rr\^2 ^t+l 2 lim 

^ \Z Kim 

avec: 
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(50) 



t+1^^ 



: Ln t+1 t+1 

Zora ~Bom 



Lp m = Ln 



Zp^ +(1-^2) Bp 



t+r 

my 



P=0,G 



Din, 



t+1 



Zqj 



t+1 



t+1 t+1 t+1 t+l t+1 

t+1 — t+1 +7 rr t+1 2 — " r- t+1 2 

_ Bom Bom 2^/2Bom 2"\/2BGm J 



t+1 



t+1 
D2m = 



r- 2 t+1 2 

V2R 



■ t+1 t+1 

Lgiii t+1 Lom t+1 
t+1 aiGm - t+1 aio m 
Bom Bom 



Dans les equations (50), on reconnaitra les parametres de phases adimensionnels B et A 
definis precedemment par les equations (19) et (20). 

Dans I'equation (49), le terme litl represente la contribution des interactions binaires entre 
constituants et ce terme est nul lorsque les coefficients d* interactions binaires de la 
representation detaillee sont tous nuls. Dans un tel cas, les coefficients d'dquilibre peuvent 

etre calcules directement a partir des parametres (EnPm) issus de I'equation d'etat 
EOS_PRO de la simulation dynamique utilisant la representation groupee. 

«t+i 

Dans le cas contraire, une expression rigoureuse du terme lim est: 

t+1 Nrb t+I Nrb ____ ____ 

(51) I^m =-rT?r X [Vai(Tm)Vaj(Tm)5ij x/^' ^ [ Vai(Tm)Vaj(T'^') 5ij y/;' ] 

Ce terme fait done intervenir les fractions molaires, dans la representation detaillee, des 

constituants dans les phases hydrocarbure au pas de temps (t+1) Xim^ Tim^ non encore 

estimees. Pour eviter une resolution implicite lourde, on peut, dans I'equation (51), 

t+1 t+I 

substituer les fractions molaires Xim, yim par les fractions molaires au pas de temps (t) 
issues des equations (46). 

C'est avec une telle approximation, non limitative de realisation, que les resultats de 
delumping : 
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- illustres Fig; 8.1 a 8.16, Fig. 9.1 a 9.16 et recapitules Fig. 10 pour une operation de 
depressurisation du fluide du cas SPE3, 

- illustres Fig. 11.1 a 11.16 et Fig. 12.1 k 12.16 et recapitules Fig.l3 pour une operation 
de vaporisation a une pression de 169 bar d'un condensat Qe condensat obtenu & 169 



ont ete obtenus. 

Dans les Fig. 8.1 a 8.16 et 9.ra 9.16 d'une part. Fig. 11.1 a 11.16 et 12.1 a 12.16 d'autre 



molaire de la phase, et, dans chaque figure, I'echelle est dilatee pour etre adaptee a la 
10 variation de la teneur du constituant dans la phase. 

L'abscisse sur les Fig. 11.1 a 11.16 et Fig. 12.1 a 12.16 est la fraction molaire de vapeur au 
cours de Foperation de vaporisation, qui croit de 0 (a I'etat initial, seule la phase 
condensat, ou huile, est presente) a 1, lorsque tout le condensat est vaporise. 

La Fig. 10 et la Fig. 13 montrent que Terreur absolue maximale obtenue sur la terieur en 
15 . chaque constituant (toujours exprimee en pourcentage molaire) reste en moyenne 
inferieure a 0.1% durant les operations de depressurisation et de vaporisation. - 

On constate done que les compositions issues de Toperation de delumping sont en bonne 
coherence avec les compositions detaillees de reference des fluides, bien que 
r approximation utilisee ne soit pas parmi les plus sophistiquees. 

20 Une autre approximation possible, toujours non limitative de realisation, est d'estimer les 

t+i t+i ^ . . . ^ . . t 

iractions moiaires Xj m , yi m a partir des coefficients d'equilibre Kim calculables a partir des 

fractions moiaires donnees par les equations (46) par: 



5 



bars lors de Foperation de depressurisation precedente). 



part, la teneur en chaque constituant dans la phase consideree est exprimee en pourcentage 



(52) 





i + (k4-i) 



m 



t+1 
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En general, les coefficients d'interactions binaires sont petits devant 1, et un calcul iteratif, 
par example par substitutions successives, ne s'avere pas necessaire, mais reste possible. 

pp) On detemiine, a chaque pas de temps (t+1), la fraction vaporisee (0^^ ) dans 
chaque maille (m), soit a partir des resultats de la simulation dynamique avec la 
representation groupee soit, si une meilleure precision est recherchee, en resolvant 
r equation de Rachford-Rice a partir des fractions molaires de chaque composant i dans la 
composition detaillee globale (Zim ) du fluide hydrocarbure dans la maille (m) au pas de 
temps (t+1) 

(53) S , = 0 



iti 1+ (Kin, -1)9 



in 



10 Les nombres de moles par phase dans chaque maille k chaque pas de temps (t+1) sont 
donnes par: 

t+1 t+1 t+1 t+1 ^ t+1 ^t+i 

(54) Non. =Nn. (l-Bn,) ; Ngi„,=N„,en. 

qq) On estime la composition detaillee de chaque phase hydrocarbure, h chaque pas de 
temps (t+1) et dans chaque maille (m), en utilisant les relations suivantes : 

t+1 t+1 

15 (55) Xim = t+r — t+i pour la phase huile 

1 + (Kim -1) 6 m 

t+1 t+1 t+1 

(56) yim — i-^l ^ . _ t+1 pour la phase gaz 

1 (Kim ~1) y m 

rr) Enfin, on obtient les nombres de moles de chaque constituant i dans la 

t+1 

representation detaillee dans chaque phase au pas de temps t+1, respectivement (Noim) 
pour la phase huile et (Ng^m) par les equations (46) remaniees en: 

i+i t+1 t+1 t+1 t+1 t+1 

20 (57) NOim = No m Xim ; Ngim = Ng m yim 

ss) Les paragraphes nn) a it) detaillent le contenu des operations pour passer du pas de 
temps t au pas de temps t+1, en supposant connus les nombres de moles de chaque 
constituant i dans la representation detaillee dans chaque phase au pas de temps t dans 
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chaque maille de la simulation avec la representation groupee, respectivement (Noim) pour 
la phase huile et (Ngim) pour la phase gaz. Pour demarrer les calculs, il suffit done de 
connaitre ces quantites a un instant ou pas de temps donne, qui peut etre en particulier 
Tinstant initial. Ceci peut etre obtenu par divers moyens connus en soi, un moyen 
consistant en un calcul d'initialisation d'une simulation compositionnelle (calcul de I'^tat 
initial uniquement) utilisant la representation compositionnelle detaillee, les entrees autres 
que compositionnelles etant les mSmes que celle de la simulation compositionnelle avec la 
representation groupee. 

Des lors, on sait decrire revolution de la composition detaillee dans chaque maille au 
cours du procede d'exploitation modelise dans un modele compositionnel utilisant une 
representation groupee grace a une memoire des parametres de Tequation d'etat de la 
representation detaillee conservee dans la methode de calcul des proprietes de phases. 

Certaines etapes de calcul peuvent apparaftre non necessaires car redondantes. Certaines^ 
de ces redondances ont Tavantage de permettre de diminuer I'espace memoire necessaire 
au stockage des donnees, mais peuvent a Tevidence etre evitees si les considerations 
d'espace memoire ne sont pas prioritaires. 

Dans le cas oii TEOS^PRO se revele valide pour le calcul des coefficients d'equilibre 
(Ki,m et t+1), le delumping des resultats d'une simulation de reservoir effectue.avec la 
methode selon I'invention, peut s'effectuer par application de la methode decrite dans le 
brevet WO/42.937 precite. 

Remarques concernant la generalisation de Vapproche 

tt) Les details foumis dans la description de la methode utilisent un formalisme 
molaire ; on aurait pu tout aussi bien utiliser un formalisme massique. 

uu) On peut remarquer que le raisonnement base sur I'hypothese maitresse «M» 
decrite dans le paragraphe d), applique a Tensemble des constituants de la representation 
de base, ensemble ci-apres designe par «B », dans les paragraphes precedents, aurait pu 
altemativement etre applique a un sous-ensemble des constituants de la representation de 
base, ci-apres designe par «R». L'ensemble complementaire.de «R» dans «B» etant 
designe par « S », les Equations (5) seraient alors reecrites (avec P=0,G) comme: 
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(58) 



ieS 



e e t e 

MMpG = ypv MMv + ypi MMi + X Ypi MMi 



MMpo = xpi MMi + MMh + E xpi MMi ; Xpi + xph = 1 - X x^ 



ypv + ypi = 1 - E ypi 



Une telle demarche peut etre interessante lorsque, pour le cas traite, il est besoin de garder 
explicites certains constituants de base dans la representation groupee. Ainsi, dans un cas 
pratique ou les conditions de reservoir ne permettraient pas de negliger la dissolution de 
certains constituants (tels le dioxyde de carbone) dans Teau, ces constituants pourtaient 
former Tensemble « S ». 

Les simples changements de variables indiques ci-apres permettent de ramener les 
equations (5) a la situation ou Ton a considere Tensemble des constituants de la 
representation de base: 



(59) 



MMpo - Z xpi MMi 



MMi 



PO 



ieS 



1 ~ S Xpi 



Xpi|^ = 



XPI 



1 - Xxpi 



XpH|j,= 



XpH 



ieS e 
1 - X Xpi 



ieS 



MMp 



MMpG-S yPiMMi 



i€S e 



, ypv ^ Ypi 

1 - Z ypi 1 - E ypi 



Les equations (6) etant resolues avec les nouvelles variables, une transformation inverse de 
10 variables permet alors d'enchainer les calculs des paragraphes g) et suivants. 

La representation groupee comporte alors les constituants (V), (I), (H) et les constituants 
de base de Tensemble « S ». Dans chaque phase de I'ensemble «E», les fractions 
molaires des constituants de base formant Tensemble «S» sont les memes dans la 
representation groupee que dans la representation de base. 

15 vv) Les compositions des phases de I'ensemble «E» dans la representation groupee 
etant ainsi obtenues, on peut suivre les grandes lignes des etapes decrites precedemment 
pour modeliser les proprietes de phases et les equilibres thermodynamiques. Des lors 
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qii'une equation" d'etat est recherchee pour la modelisation avec la representation groupee 
des proprietes de phases et/ou des equilibres thennodynamiques, Tindication est de 
conserver pour les constituants de base formant ] 'ensemble «S » les parametres de ces 
constituants dans Tequation d'etat de la description detaillee. Les seuls parametres a 
5 estimer restent alors les parametres des constituants (V), (I), (H) dans I'equation d'etat, et 
le cas echeant les parametres d'interactions binaires de ces constituants avec les 
constituants de Tensemble « S ». Lorsque ces demiers sont poses comme nuls, on peut 
suivre la ligne directrice suivie dans I'ensemble des paragraphes n) h w) qui est de 
rechercher les parainetres d'equation d'etat EOS^PRO utiles aux calculs, dans les 
10 conditions de reservoir, de densites des phases (via les facteurs de compressibilite) et de 
viscosites par des resolutions simples (non iteratives) de systemes d'equations lineaires. 



f 



1 er depot 



50 

REVENDICATIONS 

1) Methode de pseudoisation pour estimer les proprietes ou le compoitement des phases 
hydrocarbure liquide et/ou vapeur a partir de donn6es relatives a un ensemble de reference 

constitue de melanges hydrocarbones dans une serie d'etats thermodynamiques decoulant 
5 de conditions d' exploitation definies d'un . gisement souterrain d'hydrocarbures, 
caracterisee en ce que : 

- on regroupe chacun des dits melanges hydrocarbones en au moins trois constituants 
(V, I, H), aucun de ces constituants ne correspondant a une selection particuliere de 
composants ou pseudo-composants de base qui serviraient a une description 

10 corapositionnelle detaillee des fluides, en considerant que les phases gaz issues de la 

separation en conditions de surface de chacun des melanges hydrocarbones sont des 
melanges dont le troisieme constituant (H) est exclu, et que les phases huile issues de 
la separation en conditions de surface de chacun des melanges hydrocarbones sont des 
melanges dont le premier constituant (V) est exclu ; 

15 - on determine par bilan de matiere les compositions des produits de separation 
comprenant, pour les produits gazeux, au moins le premier et le deuxieme constituant 
(V, I) en proportions variables, et pour les produits liquides, au moins le deuxieme et le 
troisieme constituant (I, H) en proportions variables ; et 

- on determine la composition a au moins trois constituants de chaque melange 
20 hydrocarbone de Tensemble de reference par combinaison des produits de sa 

separation au prorata des quantit^s de chaque produit de separation. 

2) Methode selon la revendication 1, dans laquelle on regroupe chacun des melanges 
hydrocarbones en seulement trois constituants (V, I, H), les phases gaz issues de la dite 
separation etant des melanges en proportions variables du premier constituant (V) et du 

25 deuxi&me constituant (I), les phases huile issues de la dite separation etant des melanges en 
proportions variables du deuxieme constituant (I) et du troisidme constituant (H), et on 
determine la composition a trois constituants. 

3) Methode selon la revendication 1 ou 2, dans laquelle les conditions de surface sont 
celles rencontrees ou prevues lors de I'exploitation du gisement. 



1 er depot 

51 



4) Methode selbn la revendication 1 ou 2, dans laquelle les conditions de surface sont 
differentes de celles rencontrees ou prevues lors de 1' exploitation du gisement 

5) Methode selon la revendication 1 ou 2, dans laquelle le bilan de matiere est un bilan 
naassique et dans laquelle on attribue a chacun des trois constituants (V, I, H) une masse 

5 molaire apres une analyse quantitative des masses molaires des produits de separation de 
Tensemble de reference. 

6) Methode selon la revendication 1 ou 2, dans laquelle on definit les donnees necessaires 
aux calculs d'equilibre et a la modelisation des proprietes de phases dans la representation 
groupee, en utilisant les compositions des phases dans la representation group6e et des 

10 donnees connues, ou estimees, a priori portant au moins sur la densite et la viscosite de 
phases huile et gaz en equilibre appartenant a T ensemble de reference. 

7) Methode selon la revendication 6, dans laquelle, lorsque les dites donnees incluent des 
donnees compositionnelles d6taill6es des phases prealablement representees par une 
equation d'etat, on definit les parametres d'une premiere equation d'etat de la 

15 representation groupee, utile a la modelisation des proprietes de phases, en utilisant ces 
donnees compositionnelles. ^ . 

8) Methode selon la revendication 6, dans laquelle on ajuste les parametres d'une seconde 
equation d'etat de la representation groupee, utile aux calculs d'equilibre, pour reproduire 
les coefficients d'equilibre de la representation groupee. 

20 9) Methode selon la revendication 8, dans laquelle on utilise, dans la procedure 
d'ajustement pour definir les parametres de la seconde equation d*etat de la representation 
groupee utile aux calculs d'equilibre, les parametres par constituant de la representation 
groupee dans la premiere equation d'etat. • ^ 

10) Methode selon Tune des revendications precedentes, dans laquelle on determine les 
25 coefficients d'equilibre des fluides dans une representation compositionnelle detaillee, a 
partir de variables et/ou parametres intervenant dans le calcul des proprietes de phases, des 
lors que les parametres utiles aux calculs des proprietes de phases dans la representation 
groupee ont ete estimes de faijon a reproduire les parametres des phases dans I'equation 
d'etat de la description compositionnelle detaillee. 
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11) Methode selon Tune des revendications precedentes, comportant un eclatement ou 
delumping pour predire en fonction du temps et dans au moins une zone 
thermodynamique, une composition detaillee d'un fluide contenu dans un gisement 
d'hydrocarbures ou produit par au moins un puits. 

5 12) Methode selon la revendication 11 dans laquelle : 

- on repr6sente le reservoir sous la forme d'un reseau de mailles (m) dont chacune 
constitue un volume elementaire rempli de fluides sous forme de une ou plusieurs 
phases, avec au moins une phase non aqueuse ; 

- on definit, pour cheque zone ou domaine thermodynamique, les fluides par une 
10 representation detaillee de base, de fafon a determiner la quantite de chaque 

constituant de base (i) dans chaque phase hydrocarbure dans chaque maille (m) a 
r instant defini comme initial pour le calcul de delumping ; 

- par zone thermodynamique pour laquelle on adopte une representation groupee des 
fluides, on determine une equation d'etat construite, prealablement a la simulation 

15 dynamique de reservoir avec la representation groupee, pour reproduire les parametres 

de phases, dans Fequation d'etat de la representation detaillee, des fluides 
hydrocarbones au cours de chemins thermodynamiques tenus pour representatifs de 
ceux que vont suivre les fluides hydrocarbones pendant la simulation maillee ; 

- on realise a un pas de temps t une simulation compositionnelle avec un nombre limite 
20 de constituants ou les propriet^s des phases sont calculees par une equation d'etat, 

ladite simulation permettant de calculer au moins dans chaque maille (m) et a des pas 
de temps consecutifs, une pression pour une phase hydrocarbure, la temperature 
lorsque celle-ci varie, les debits des phases entre mailles et au niveau des perforations 
de production et d'injection, et les valeurs de parametres et/ou proprietes de phase 
25 intervenant dans F expression formelle des coefficients d'equilibre de la representation 

detaillee, et on memorise ces diverses quantites ; 

- on estime au pas de temps suivant (t+1) la fraction molaire de chaque constituant i 
dans la composition detaillee globale du fluide hydrocarbone dans la maille (m) par 
bilan matiere sur la maille (m) ; 



ler depot 



,53 

on determine en utilisant les quantites memorisees, au mSme pas de temps (t+1) et 
dans chaque inaille (m), les coefficients d'equilibre de chaque constituant (i) dans la 
representation detaillee ; 

on determine, au meme pas de temps (t+1), la fraction vaporisee dans chaque maille 
(m);et 

on estime la composition detaillee de chaque phase hydrocarbonee, au meme pas de 
temps (t+1) et dans chaque maille (m). 

13) Methode selon Tune des revendications precedentes, ou differentes etapes sont 
traduites de fa9on h, produire des resultats intermediaires ou finaux utilisables dans des 
equations exprimees selon un formalisme molaire ou massique. 
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FIG.7A 
Fraction molaire de C1 dans le gaz 
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FIG.7-1 
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FIG.7-4 
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FIG.7-7 

Fraction molaire de C6 dans le gaz 
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FIG.7-10 
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